PHOTOGRAPH  THIS  SHEET 


%■ 


QC 


LEVEL  INVElwRY 

6ENCRAL  ELECTRIC  CO  CINCINNATI  OHZ9  AIRCRAFT 

ENA  INC  GROUP?  ncrr  i  ng  tTIOP  AMP 

high  velocity  jet  NOISE*  SOURCE  LOCATION  ** 

REDUCTION^  task  3.  experimental  investigation  of  suppression 

PRINCIPLES,  VOLUME  IV  -  LASER  VELOCIMETER  TIME  DEPENDENT  CROSS 
CORRELATION  MEASUREMENTS,  FINAL  REPT,  DEC, '78, 


H  REPT,  NO,  R78AEG627  CONTRACT  DOT-0S-30034  FAA.RD-76-79-3-4, 


tribution  statement  a 


Approved  foi  public  releosai 
_ Distribution  Unlimited 


DISTRIBUTION  STATEMENT 


ACCESSION  FOR 


NTIS  GRAAI 

DTIC  TAB 

UNANNOUNCED 
JUSTIFICATION 


BV 


DISTRIBUTION  / 


AVAILABILITY  CODES 


BIST  AVAIL  AND/OR  SPECIAL 


J 

s 


DTIC 

ELECTE 
JAN  29  1981 


DATE  ACCESSIONED 


DISTRIBUTION  STAMP 


81  1  27  002 


DATE  RECEIVED  IN  DTIC 

PHOTOGRAPH  THIS  SHEET  AND  RETURN  TO  DTIC-DDA-2 


.  *  FORM  . 

DTIC  OCT  79  70A 


DOCUMENT  PROCESSING  SHEET 


REPORT  NO.  FAA-RD-76-79,  III  -  IV 


R78AEG627 


HIGH  VELOCITY  JET  NOISE 
SOURCE  LOCATION  AND  REDUCTION 

TASK  3  -  EXPERIMENTAL  INVESTIGATION 
OF  SUPPRESSION  PRINCIPLES 

Volume  IV  •  Laser  Velocimeter  Time  Dependent 
Cross  Correlation  Measurements 

AUTHORS : 

P.R.  Knott 
P.F.  Scott 
P.W.  Mossey 

GENERAL  ELECTRIC  COMPANY 
AIRCRAFT  ENGINE  GROUP 
CINCINNATI,  OHIO  4*215 


DECEMBER  1978 
FINAL  REPORT 


Document  is  available  to  the  U.S.  public  through 
the  National  Technical  Information  Service, 
Springfield,  Virginia  22161 . 


Prepared  for 

U.S.  DEPARTMENT  OF  TRANSPORTATION 

FEDERAL  AVIATION  ADMINISTRATION 
Systems  Research  &  Development  Service 
Washington,  D.C.  20590 


1.  Report  No.  2.  Government  Accession  No. 

FAA-RD-76-79,  III  -  IV 


4.  Title  and  Subtitle 

High  Velocity  Jet  Noise  Source  Locution  and  Reduction  Task  3  - 

Experimental  Investigation  of  Suppression  Principles 

Volume  IV  -  Laser  Velocimeter  Time  Dependent  Cross  Correlation 


7.  Author(s) 
P.R.  Knott 
P.F.  Scott 


P.W.  Mossey 


9.  Performing  Organization  Name  and  Address 

General  Electric  Company 
Group  Advanced  Engineering  Division 
Aircraft  Engine  Group 
Cincinnati,  Ohio  4521S 


12.  Sponsoring  Agency  Name  end  Address 

U.S.  Department  of  Transportation, 
Federal  Aviation  Administration 
Systems  Research  and  Development  Service 


3.  Recipient's  Catalog  No. 


5.  Report  Data 
December  1978 


6.  Per'otming  Organization  Code 


8.  Performing  Organization  Report  No. 
R78AEC627 


10.  Work  Unit  No. 


11.  Contract  or  Grant  No. 
D0T-0S-30034 


13.  Type  of  Report  and  Period  Covered 

Technical  Report 
October  1974-October  1977 


14.  Sponsoring  Agency  Coda 
ARD-S50 


15.  Supplementary  Notes  g,s,  Zuckerman,  Program  Manager;  this  report  is  in  partial  fulfillment  of  the 

subject  program.  Related  documents  to  be  issued  in  the  course  of  the  program  include  final  reports 
of  the  following  .Tasks :  Task  1  -  Activation  of  Facilities  and  Validation  of  Source  Location 
Techniques;  Task  1  Supplement  -  Certification  of  the  General  Electric  Jet  Noise  Anechoic  Test 
Facility;  Task  2  -  Theoretical  Developments  and  Basic  Experiments;  Task  4  -  Development/Evaluation 
of  Techniques  for  "Inflight"  Investigation;  Task  5  -  Investigation  of  "Inflight"  Aero-Acoustic 
Effects;  Task  6  -  Noise  Abatement  Nozzle  Design  Guide. 


16.  Abstract  Experimental  investigations  were  conducted  of  suppression  principles;  including  developing 
an  experimental  data  base,  developing  a  better  understanding  of  jet  noise  suppression  principles, 
and  formulating  empirical  methods  for  the  acoustic  design  of  jet  noise  suppressors.  Acoustic 
scaling  has  been  experimentally  demonstrated,  and  five  "optimum"  nozzles  have  been  selected  for 
subsequent  anechoic  free-jet  testing. 

In-jet/in-jet  and  in- jet/far-field  exhaust  noise  diagnostic  measurements  were  conducted  using  a 
Laser  Velocimeter  (LV).  Measurements  were  performed  on  a  conical  nozzle  ai:d  a  coannular  plug 
nozzle.  Two-point,  space/time  measurements  using  a  two-LV  system  were  completed  for  the  conical 
nozzle.  Measurements  of  mean  velocity,  turbulent  velocity,  eddy  convection  speed,  and  turbulent 
length  scale  were  made  fer  a  subsonic  ambient  jet  and  for  a  sonic  heated  jet.  For  the  coannular 
plug  nozzle,  a  similar  series  of  two-point,  laser-correlation  measurements  were  performed.  In 
addition,  cross  correlations  between  the  laser  axial  component  of  turbulence  and  a  far-field 
acoustic  microphone  were  performed. 

This  volume  is  part  of  the  four  volume  set  that  constitutes  the  Task  3  final  report.  The  other 
volumes  are: 

Volume  I  -  Verification  of  Suppression  Principles  and  Development  of  Suppression 
Prediction  Methods 

Volume  II  -  Parametric  Testing  and  Source  Measurements 

Volume  III  -  Suppressor  Concepts  Optimization 


17.  Key  Words  (Suggested  by  Author(i|  I 

Laser  Velocimeter,  Real-Time,  Velocity 
Measurements,  Cross-Correlation 


18.  Distribution  Statement 

Document  is  available  to  the  public  through 
the  National  Technical  Information  Center 
Springfield,  Virginia 


19.  Security  Gessif.  (of  this  report) 

20.  Security  Claud.  (ol  this  page) 

21.  No.  of  Pages 

Unclassified 

Unclassified 

3  36 

*For  sale  by  the  National  Technical  Information  Service,  Springfield,  Virginia  22151 


NASAC-161  *•?!) 


METRIC  CONVERSION  FACTORS 


PREFACE 


This  report  describes  the  work  performed  under  Task  3  of  the  DOT/FA A 
High  Velocity  Jet  Noise  Source  Location  and  Reduction  Program  (Contract  DOT- 
OS-  30034).  The  objectives  of  the  contract  were: 

•  Investigation,  including  scaling  effects,  of  the  aerodynamic  and 
acoustic  mechanisms  of  various  jet  noise  suppressors. 

•  Analytical  and  experimental  studies  of  the  acoustic  source  distri¬ 
bution  in  such  suppressors,  including  identification  of  source 
location,  nature,  and  strength  and  noise  reduction  potential. 

•  Investigation  of  in-flight  effects  on  the  aerodynamic  and  acoustic 
performance  of  these  suppressors. 

The  results  of  these  investigations  are  expected  to  lead  to  the  preparation 
of  a  design  report  for  predicting  the  overall  characteristics  of  suppressor  con¬ 
cepts,  from  models  to  full  scale,  static  to  in-flight  conditions,  as  well  as 
a  quantitative  and  qualitative  prediction  of  the  phenomena  involved. 

The  work  effort  in  this  program  was  organized  under  the  following  major 
Tasks,  each  of  which  is  reported  in  a  separate  Final  Report: 

Task  1  -  Activation  of  Facilities  and  Validation  of  Source  Location 
Techniques. 

Task  2  -  Theoretical  Developments  and  Basic  Experiments. 

Task  3  -  Experimental  Investigation  of  Suppression  Principles. 

Task  4  -  Development  and  Evaluation  of  Techniques  for  "In-flight" 
Investigation. 

Task  5  -  Investigation  of  "In-flight"  Aero-Acoustic  Effects  on  Suppressed 
Exhausts. 

Task  6  -  Preparation  of  Noise  Abatement  Nozzle  Design  Guide  Report. 

Task  1  was  an  investigative  and  survey  effort  designed  to  identify 
acoustic  facilities  and  test  methods  best  suited  to  jet  noise  studies. 

Task  2  was  a  theoretical  effort  complemented  by  theory  verification  experi¬ 
ments  which  extended  across  the  entire  contract  period  of  performance. 

The  subject  of  the  present,  Task  3  report  series  (FAA-RD-76-79III  —  I, 

II,  III,  and  IV)  was  formulated  as  a  substantial  part  of  the  contract 
effort  to  gather  various  test  data  on  a  wide  range  of  high  velocity  jet 
nozzle  suppressors.  These  data,  together  with  supporting  theoretical  advances 
from  Task  2,  have  led  to  a  better  understanding  of  jet  noise  and  jet  noise 


lil 


suppression  mechanisms,  as  well  as  to  a  validation  of  scaling  methods.  Task 
3  helped  to  identify  several  "optimum"  nozzles  for  simulated  in-flight 
testing  under  Task  5,  and  to  provide  an  extensive,  high  quality  data  bank 
leading  to  formulation  of  methods  and  techniques  useful  for  designing  jet 
noise  suppressors  for  application  in  the  Task  6  design  guide  as  well  as  in 
future  studies. 

Task  4  was  similar  to  Task  1,  except  that  it  dealt  with  the  specific 
test  facility  requirements,  measurement  techniques,  and  analytical  methods 
necessary  to  evaluate  the  "in-flight"  noise  characteristics  of  simple  and 
complex  suppressor  nozzles.  This  effort  provided  the  capability  to  conduct 
the  "flight"  effects  test  program  of  Task  5. 
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Velocity  component  perpendicular  to  the  fringes 

Distance  along  cptic  axis  of  lens  (Sections  4  and  5);  coordinate 
along  jet  axis 

Position  vector  from  origin  to  observer 
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Position  vector  from  origin  to  jet  noise  source  element 

The  point  on  the  unit  Gausian  distribution  such  that  the  area 
under  the  distribution  is  a: 
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Included  angle  between  laser  beams;  also  see  definition  of  Za 
Slit  half  angle  of  subtention 
Unity  sensor 

Optical  path  length  (Sections  4  and  5);  the  width  of  the  time  grid 
used  to  reconstruct  the  correlation  function  (Sections  7  and  8) 

Scattering  frequency 

Observation  angle 

Convection  amplification  factor  (see  Equation  6-5) 

Wavelength  of  coherent  radiation  (Sections  4  and  5);  the  average 
rate  of  arrival  of  particles  at  the  probe  volume  of  the  laser 
ve  locimeter 

The  unit  step  function 

U(t)  =  l  ;  t  >_  o 
u(t)  =  0  ;  t  <  o 
The  Dirac  delta  function 

The  expected  value  or  mean  of  the  random  variable  x. 

The  operation  distance  between  LV  probe  volumes  in  a  two-point 
velocity  correlation  measurement 

Fluid  density 

The  standard  deviation  of  the  random  variable  x;  also  the  turbu¬ 
lence  level 

Time 
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The  retarded  time  (r/ac)  J 

The  value  of  t  where  |r  (t) I  is  maximum  ! 
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The  viscous  stress  tensor 


Typical  radian  frequency  of  turbulence 
Radian  frequency  in  the  turbulence 
Weight  flow 

Scripts 

Inner  stream  of  coannular  plug  nozzle 
Outer  stream  of  coannular  plug  nozzle 


1.0  SUMMARY 


The  High  Velocity  Jet  Noise  Source  Location  and  Reduction  Program  (Con¬ 
tract  DOT-OS-30034)  was  conceived  to  bring  analytical  and  experimental  know¬ 
ledge  to  bear  on  understanding  the  fundamentals  of  jet  noise  for  simple  and 
complex  suppressors. 

Task  3,  the  subject  of  this  report,  involved  the  experimental  investi¬ 
gation  of  suppression  principles,  including  developing  an  experimental  data 
base,  developing  a  better  understanding  of  jet  noise  suppression  principles, 
and  formulating  empirical  methods  for  the  acoustic  design  of  jet  noise  sup¬ 
pressors.  Acoustic  scaling  has  been  experimentally  demonstrated,  and  five 
"optimum"  nozzles  were  selected  for  anechoic,  free-jet  testing  in  Task  5. 

Volume  I  -  Verification  of  Suppression  Principles  and  Development  of 
Suppression  Prediction  Methods  -  Some  of  the  experimental  studies  (reported 
in  Volume  II)  involved  acquisition  of  detailed,  far-field,  acoustic  data  and 
of  aerodynamic  jet-flow-field  data  on  several  baseline  and  noise-abatement 
nozzles.  These  data  were  analyzed  and  used  to  validate  the  theoretical  jet 
noise  prediction  method  of  Task  2  (referred  to  as  M*G*B,  designating  the 
authors'  initials)  and  to  develop  and  validate  the  empirical  noise-prediction 
method  presented  herein  (referred  to  as  M*S,  designating  the  last  name  ini¬ 
tials  of  the  authors).* 

The  Task  2  theoretical  studies  conclude  that  four  primary  mechanisms  in¬ 
fluence  jet  noise  suppression:  fluid  shielding,  convective  amplification, 
turbulent  mixing,  and  shock  noise.  A  series  of  seven  suppressor  configura¬ 
tions  (ranging  from  geometrically  simple  to  complex)  were  evaluated  in  Task  3 
to  establish  the  relative  importance  of  each  of  the  four  mechanisms.  Typical 
results  of  this  evaluation  of  noise  mechanisms  are  summarized  in  Figure  1-1 
in  terms  of  perceived  noise  level  ( PNL)  directivity  for  a  conical  nozzle.  In 
general,  mechanical  suppressors  exhibit  a  significant  reduction  in  shock 

*The  Task  3  empirical  (M*S)  method  was  initially  intended  for  nozzle 
geometries  which  could  not  be  modeled  in  the  purely  analytical  Task  2 
(M*G*B)  method  (a  multielement  nozzle  with  a  treated  ejector,  for 
example ) . 
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noise  relative  to  a  baseline  conical  nozzle,  reduce  the  effectiveness  of 
fluid  shielding  (increase  rather  than  suppress  noise),  reduce  the  effective¬ 
ness  of  convective  amplification  (reduce  noise),  and  produce  a  modest  reduc¬ 
tion  in  turbulent  mixing  noise.  The  largest  amount  of  shock  noise  reduction 
correlates  with  the  suppressor  which  has  the  smallest  characteristic  dimension. 
Fluid  shielding  decreases  because  suppressors  cause  the  mean  velocity  and 
temperature  of  the  jet  plume  to  decay  faster  than  the  conical  baseline.  A 
reduction  in  convection  Mach  number  (and  hence  in  convective  amplification) 
occurs  because  a  suppressor  plume  decays  very  rapidly.  Turbulent  mixing 
noise  is  reduced  through  alteration  of  the  mixing  process  that  results  from 
segmenting  the  exhaust  jet. 

Aerodynamic  flow-field  measurements  (mean-velocity  profiles)  were  demon¬ 
strated  to  be  useful  in  verifying  the  flow-field  predictions  which  were  cal¬ 
culated  by  the  M*G*B  (theoretical)  noise-prediction  program.  Noise  source 
location  devices  such  as  the  Ellipsoidal  Mirror  (EM)  were  demonstrated  to  be 
less  useful  than  the  Laser  Velocimeter  (LV)  for  the  M*G*M  theory  verification 
studies  because  the  LV  provides  data  which  may  be  directly  compared  with 
predictions  made  using  the  M*G*M  program.  Axial  and  radial  mean-velocity 
profiles  are  typical  examples  of  such  comparisons. 

The  empirical  M*S  jet  noise  prediction  method  has  been  developed  to  pre¬ 
dict  the  static  acoustic  characteristics  of  multielement  suppressors  appli¬ 
cable  to  both  advanced  turbojets  and  variable-cycle  engines  (which  are  repre¬ 
sentative  of  power  plants  for  future  supersonic  cruise  aircraft).  The  effect 
of  external  flow  on  the  M*S  jet  noise  prediction  is  discussed  in  the  Task  6 
Design  Guide  Report.  Inputs  required  to  use  the  M*S  computational  procedure 
include:  element  type,  element  number,  suppressor  area  ratio  and  radius 
ratio,  chute-spoke  planform  and  cant  angle,  and  plug  diameter.  The  predic¬ 
tion  accuracy  is  estimated  to  be  _+3 . 3  Effective  Perceived  Noise  Decibels 
(EPNdB)  at  a  95%  confidence  level.  Figure  1-2  illustrates  the  correlation 
between  measured  and  predicted  EPNLs  for  all  types  of  suppressors. 

The  merits  of  both  the  M*S  and  M*G*B  computational  techniques  can  be 
stated  as  follows.  The  empirical  (M*S)  jet  noise  prediction  method,  based  on 
correlations  of  scale-model  jet  data,  serves  as  a  useful  preliminary  design 


3 


Measured  EPNL.  Relative  to  Reference,  M,  EPNdB 


•  Flyover  calculation  using  static  data  corrected  to  free-field  conditions. 


•  The  Reference"  level  is  the  predicted  value  of  noise  for  each  nozzle, 
at  a  specified  set  of  thermodynamic  conditions,  plus  an  arbitrary 
value  of  100  dB. 


Figure  1-2.  Correlation  Between  Measured  and  Predicted  Effective  Perceived 
Noise  Level,  EPNL,  for  all  Types  of  Suppressor  Nozzles. 


and  prediction  tool  for  selecting  the  basic  nozzle  type  (chute,  spoke,  multi¬ 
tube,  etc.)  and  primary  geometric  parameters  (element  number,  area  ratio, 
etc.)  for  a  given  application.  It  is  also  useful  in  evaluating  the  acoustic 
performance  of  a  given  suppressor  nozzle,  provided  the  nozzle  is  one  of  the 
types  from  which  the  correlation  was  derived.  Further,  the  method  is  useful 
for  doing  parametric  studies  since  the  computation  procedure  is  relatively 
simple  and  economical  of  both  computer  time  and  cost.  The  theoretical 
(M*G*B)  prediction  method,  on  the  other  hand,  is  more  suited  to  detailed  de¬ 
sign  and  analysis  of  a  suppressor  nozzle.  It  can  supply  detailed  information 
on  the  jet  plume  flow  development  as  well  as  the  far-field  acoustic  character¬ 
istics.  It  is  also  capable  of  evaluating  changes  in  nozzle  planform  shape, 
element  placement  and  spacing,  etc.  In  addition,  the  theoretical  prediction 
model  is  a  useful  diagnostic  tool,  capable  of  assessing  the  relative  roles  the 
various  mechanisms  play  in  the  noise  suppression  process,  and  can  also  serve 
as  a  source  location  analysis  tool. 

Volume  II  -  Parametric  Testing  and  Source  Measurements  -  A  parametric 
experimental  series  was  conducted  to  provide  far-field  acoustic  data  on  47 
baseline  and  suppressor  nozzle  conf igurat ions  and  to  provide  aerodynamic 
nozzle  performance  on  18  of  the  configurations.  The  data  presented  in  this 
volume  were  taken  for  use  in  the  current  program  as  well  as  to  provide  an  ex¬ 
tensive,  high-quality,  data  base  for  future  studies.  The  impact  of  varying 
the  area  ratio  and  velocity  ratio  of  dual-flow,  baseline  nozzle  configurations 
was  investigated,  and  the  importance  of  shock  noise  was  assessed.  The  impact 
of  varying  area  ratio  and  element  number  was  parametrically  studied  for  both 
single  and  dual-flow  suppressors;  core  plug  geometry,  velocity  ratio,  and 
weight  flow  ratio  were  evaluated  for  dual-flow  suppressors.  These  studies 
establish  absolute  static  suppression  levels  on  the  basis  of  normalized  maxi¬ 
mum  PNL,  for  several  families  of  suppressor  nozzles,  as  illustrated  in  Fig¬ 
ure  1-3. 

Parametric  testing  identified  the  following  primary  trends  for  single¬ 
flow  and  for  dual-flow  suppressors  during  static  operation: 
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Figure  1-3.  Typical  Peak  Static  Noise  Suppression  Characteristics 


Single  Flow 

•  Suppression  increases  with  increasing  area  ratio  at  high  jet 
velocity. 

•  Suppression  decreases  with  increasing  area  ratio  at  low  jet 
velocity. 

•  Suppression  level  is  affected  by  element  type  (spoke  systems 
suppress  slightly  better  than  chutes). 

Dual  Flow 

•  Suppression  increases  with  increasing  area  ratio.. 

•  Suppression  increases  with  increasing  element  number  at 
high  jet  velocity. 

•  Suppression  level  is  affected  by  core  plug  geometry  [by  2 
to  3  decibels  (dB)]. 

•  Suppression  increases  3  to  4  dB  when  a  treated  ejector  is 
added  to  a  suppressor  configuration. 

Selective,  free-jet  tests  conducted  on  eight  configurations  indicate 
that  suppression  generally  decreases  in  flight.  Typical  static  versus  free- 
jet  results  are  shown  in  Table  1-1. 

The  aerodynamic  performance  test  recorded  on  18  of  the  configurations 
at  both  static  and  wind-on  conditions  is  also  included  in  this  volume.  Base 
pressure  measurements  were  taken  on  several  of  the  models  in  order  to  deter¬ 
mine  base  drag  (which  is  thought  to  be  responsible  for  the  poor  aerodynamic 
performance  of  most  mechanical  suppressors  in  flight).  These  wind  tunnel 
tests  identified  the  following  primary  trends  in  aerodynamic  performance: 

•  Performance  decreases  with  increasing  element  number. 

•  Performance  increases  with  increasing  chute  depth. 

•  Performance  increases  with  increasing  ratio  of  inner  flow  area 
to  outer  flow  area. 

•  Performance  is  affected  by  element  type  (chutes  perform  better 
than  spokes  because  spokes  have  higher  base  drag). 
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Table  1-1.  Typical  Summary  of  Nozzle  Static  and  Projected 
Flight  Peak  PNL  Suppression  Characteristics. 


•  Suppression  Levels  are  Relative  to  a  Conical 
Nozzle  at  Equivalent  Flight  Conditions 

•  Vj  =  2500  ft/sec 


Configuration 

Suppression 

Static 

Level , 

Flight 

Plug  Nozzle  -  0.789  Radius  Ratio 

1.3 

3.0 

Plug  Nozzle  -  0.85  Radius  Ratio 

2.3 

3.7 

8-Lobe  Nozzle 

5.6 

5.6 

AR  =  2.5  36-Chute  Nozzle 

13.5 

10.9 

AR  *  2.5  36-Chute  Nozzle  with  Auxiliary 
Flow 

12.5 

9.4 

104-Tube  Nozzle 


12.0 


12.0 


The  base  pressure  correlations  provide  a  procedure  for  predicting  sup¬ 
pressor  nozzle  aerodynamic  performance. 

Volume  III  -  Suppressor  Concepts  Optimization  -  Several  studies  were  con¬ 
ducted  to  attempt  an  optimization  of  suppressor  concepts.  The  end  product  of 
this  overall  effort  was  to  design  five  nozzles  for  static  and  free-jet  testing 
in  Task  5.  Trade  studies  of  performance  versus  suppression,  aircraft  inte¬ 
gration  studies,  and  development  of  a  figure  of  merit  method  of  analysis  all 
make  up  the  activities  in  this  "optimization"  process. 

Trade  studies  of  suppression  versus  aerodynamic  performance  indicate  that 
a  properly  selected  and  designed  mechanical  suppressor  can  attain  a  delta 
suppression  to  delta  thrust  coefficient  ratio  (aPWL/aCj  )  of  almost  3.0 
(based  on  static  suppression  and  wind-on  aerodynamic  performance). 

The  aircraft  integration  study  consisted  of  ranking  nine  baseline  and 
suppressor  nozzles  with  respect  to  performance  level,  suppression  level, 
weight,  impact  on  aircraft  mission  range,  and  noise  footprint.  In  general, 
suppression  level  was  found  to  be  the  most  important  design  variable,  with 
performance  and  weight  ranking  second  and  third,  respectively. 

The  appropriate  figure  of  merit,  considering  all  the  design  variables, 
was  found  to  be  aircraft  range.  However,  use  of  range  as  the  figure  of  merit 
requires  that  the  aircraft  mission  be  specified,  and  several  techniques  for 
cursorily  ranking  the  suppressors  based  solely  on  suppression  level,  perfor¬ 
mance,  and  weight  may  also  be  identified,  k  summary  of  the  range  versus  noise 
characteristics  of  typical  nozzle  configurations  is  presented  in  Figure  1-4. 
Once  a  noise  goal  is  specified,  adding  a  suppressor  provides  a  significant 
range  improvement  over  an  unsuppressed  system  because  adding  a  suppressor  is 
less  costly  than  reducing  noise  by  enlarging  the  engine  to  reduce  jet  velocity. 

The  design  of  the  five  optimum  nozzles  was  based  on  data  from  previous 
studies,  performed  by  government  and  industry,  on  the  M*G*B  and  M*S  models 
discussed  above  and  on  the  parametric  data  obtained  in  the  acoustic  and  aero¬ 
dynamic  performance  test  series  reported  in  Volume  II.  The  configurations  were 
designed  and  fabricated  for  open-throat,  anechoic,  free-jet  testing  in  Task  5. 
The  configurations  chosen  for  evaluation  were:  (1)  a  32-chute,  single-flow 
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Summary  of  Range  and  Noise  Characteristics  for  Several  Baseline  and  Suppression 
Levels. 


nozzle;  (2)  a  40-shallow-chute,  dual-flow  nozzle;  (3  and  4)  a  36-chute,  dual¬ 
flow  nozzle,  with  and  without  a  treated  ejector;  and  (5)  a  54-element,  co- 
planar-mixer,  plug  nozzle. 

Demonstration  of  acoustic  scaling  for  several  suppressor  configurations 
was  conducted  to  assure  the  adequacy  of  using  scale-model  results  to  project 
full-scale  suppression  levels.  Full-scale  data  were  obtained  on  several  sup¬ 
pressor  conf igurations  using  J79  and  J85  engines.  The  suppressors  evaluated 
were:  (1)  a  baseline  conical  nozzle,  (2)  a  32-chute  nozzle  with  and  without 
a  treated  ejector,  (3)  an  8-lobe  nozzle,  and  (4)  a  104-tube  nozzle.  Scale- 
model  data  were  obtained  for  these  same  configurations  to  allow  comparison  of 
scale-model  and  full-scale  results.  In  general,  peak  full-scale  suppression 
levels  projected  from  scale-model  data  were  verified  by  the  full-scale  engine 
results.  Directivity  patterns  were  duplicated  within  +2  PNdB  (the  largest 
differences  occurring  with  the  conical  nozzle  configuration).  Some  spectral 
anomalies  were  observed  for  select  cases;  however,  they  were  not  of  suffici¬ 
ent  magnitude  to  invalidate  the  scale-model  results.  The  conclusion  re¬ 
sulting  from  this  study  is  that  full-scale  noise  levels  can  be  predicted  from 
scale-model  test  results  using  Strouhal  scaling  laws. 

Volume  IV  -  Laser  Velocimeter  Time  Dependent  Cross  Correlation  Measure¬ 
ments  -  In-jet/in-jet  and  in-jet/f ar-f ield  exhaust  noise  diagnostic  measure¬ 
ments  conducted  using  a  Laser  Velocimeter  (LV)  are  reported  in  this  volume. 
Measurements  were  performed  on  a  conical  nozzle  and  a  coannular  plug  nozzle. 
Two-point,  space/time  measurements  using  a  two-LV  system  were  completed  for 
the  conical  nozzle.  Measurements  of  mean  velocity,  turbulent  velocity,  eddy 
convection  speed,  and  turbulent  length  scale  were  made  for  a  subsonic  ambient 
jet  and  for  a  sonic  heated  jet.  For  the  coannular  plug  nozzle,  a  similar 
series  of  two-point,  laser-correlation  measurements  were  performed.  In  addi¬ 
tion,  cross  correlations  between  the  laser  axial  component  of  turbulence  and 
a  far-field  acoustic  microphone  were  performed. 

Volumes  I,  II,  III,  and  IV  contain  the  results  of  a  comprehensive  effort 
to  identify  and  integrate  the  theoretical  studies,  parametric  test  data, 
acoustic  and  performance  diagnostic  measurements,  and  system  studies.  A 
logical  procedure  has  evolved  for  conducting  suppressor  design  trade-offs. 


2.0  INTRODUCTION 


The  first  20  years  of  commercial  aircraft  operation  with  jet  propulsion 
have  clearly  demonstrated  the  need  for  effective  high  velocity  jet  noise 
suppression  technology  in  order  to  meet  community  acceptance.  Aircraft  sys¬ 
tem  studies  show  that  an  efficient  jet  noise  suppression  device  is  required 
if  a  commercial  supersonic  aircraft  is  to  be  economically  viable  as  well  as 
environmentally  acceptable.  The  current  state  of  the  art  of  high  velocity 
jet  noise  suppression  would  make  a  supersonic  transport  (SST),  with  advanced 
technology  engines,  meet  1969  noise  rules  (at  best).  This  state  of  the  art 
is  represented  by  the  material  in  References  1  through  6. 

Reference  1  describes  analytical  and  experimental  investigations  which 
were  conducted  in  the  early  1960's.  This  study  established  a  basis  for  de¬ 
velopment  of  mathematical  and  empirical  methods  for  the  predictions  of  jet- 
flow-field,  aerodynamic  characteristics  and  for  determining  the  directional 
characteristics  of  jet  noise  suppressors.  This  work  was  limited  in  the  sense 
that  the  suppressors  evaluated  had  only  modest  suppression  potential,  and  the 
measurement  techniques  available  did  not  allow  the  acquisition  of  high- 
frequency,  spectral  data  necessary  to  establish  full-scale,  PNL  suppression 
levels . 

The  development  of  commercial  SST  vehicles  by  the  U.S.  and  by  the  British- 
French  multinational  corporation  in  the  1960's  placed  extreme  emphasis  on  the 
need  for  effective  and  efficient  noise  suppression  devices.  Phase  I  of  work, 
conducted  by  the  Boeing  and  General  Electric  companies,  is  summarized  in 
References  2  and  3.  Primary  emphasis  was  on  jet  noise  suppressor  development 
through  model  and  engine  testing  applicable  to  an  afterburning  turbojet 
engine.  Suppressor  designs  were  based  primarily  on  empirical  methods.  Phase 
II  of  this  effort,  References  4  and  5,  continued  the  suppressor  development 
with  a  stronger  emphasis  placed  on  the  integration  of  analytical  studies  and 
experimental  test  data.  Specifically,  the  Boeing  Company  concentrated  on 
optimization  of  tube-type-suppressor  systems  and  related  semiempir ical  pre¬ 
diction  methods.  General  Electric  focused  on  the  development  both  of  chute 
and  of  tube-type-suppressor  systems  with  primary  emphasis  placed  on  optimiza¬ 
tion  of  chute-type-suppressor  nozzles. 


Similar  studies  were  conducted  by  the  British  and  French  in  development 
of  the  Concorde,  and  typical  results  are  summarized  in  Reference  6. 

The  design  technology  represented  in  References  1  through  6  is  primarily 
semiempir ical .  The  absence  of  general  design  rules  based  on  engineering 
principles  led  to  the  Government's  formulation  of  the  High  Velocity  Jet  Noise 
Program,  Contract  DOT--OS-30034 ,  in  1973.  The  purpose  has  been  to  achieve 
fundamental  understanding,  on  a  quantitative  basis,  of  the  mechanisms  of  jet 
noise  generation  and  suppression  and  to  develop  design  methods. 

This  report  presents  the  results  of  Task  3  of  the  contract.  It  provides 
the  experimental  data  base  which  was  used  in  conjunction  with  the  supporting 
theories  from  Task  2  to  develop  a  better  understanding  of  jet  noise  and  jet 
noise  suppression. 

The  report  is  organized  into  four  volumes  ( FAA-RD-76-79 ,  III  -  I,  II, 

III,  IV)  and  is  presented  in  a  format  consistent  with  the  Task  3  work  plan 
division  of  the  subtasks.  Volume  I  is  entitled  "Verification  of  Suppression 
Principles  and  Development  of  Suppression  Prediction  Methods."  Volume  II  is 
a  data  report  entitled  "Parametric  Testing  and  Source  Measurements,"  and 
Volume  III  is  an  analysis  report  entitled  "Suppressor  Concepts  Optimization," 
Volume  IV,  under  this  cover,  is  an  analysis  report  entitled  "Laser  Velocimeter 
Time  Dependent  Cross  Correlation  Measurement." 

Volume  I  uses  the  data  base  (Volume  II)  and  the  Task  2  theoretical  model 
(Reference  7)  to  postulate  the  suppression  mechanisms.  Volume  I  also  presents 
an  independent,  empirical,  static  jet-noise-prediction  method  which  was  de¬ 
veloped  from  engineering  correlations  of  the  test  data.  Volume  II  presents 
the  data  and  results  of  the  parametric  acoustic  tests,  the  aerodynamic  per¬ 
formance  tests,  and  the  Laser  Velocimeter  tests.  Volume  III  presents  the 
results  of  a  trade  study  of  performance  versus  suppression,  an  aircraft  inte¬ 
gration  study,  a  "figure  of  merit"  methodology,  and  a  summary  of  the  five 
"optimum"  nozzles  selected  for  testing  in  Task  5.  An  acoustic-scaling  investi¬ 
gation  was  conducted  to  support  the  suppressor  concepts  optimization  activi¬ 
ties  and  is  presented  as  an  Appendix  to  Volume  III.  Volume  IV  presents  the 
results  of  the  in- jet/ in- jet  and  in-jet/ far-field  cross  correlation  investiga¬ 
tions  . 
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The  work  reported  in  the  present  volume  describes  a  number  of  new  laser 
velocimeter  processor  and  data  handling  concepts  that  are  required  for  real¬ 
time  cross-correlation  type  measurements.  Of  particular  note  are  discussions 
of  a  filter  bank  laser  velocimeter  processor  design  approach,  combined  with 
the  use  of  a  quantized  product -factoring  computational  concept  for  data 
handling  that  has  enabled  the  successful  measurement  of  two-point,  space-time, 
cross-correlation  measurements  in  realistic  gas  flows.  A  series  of  laser 
velocimeter  cross-correlation  measurements  (using  two  laser  velocimeters)  are 
reported  for  a  conical  nozzle  and  a  coannular  plug  nozzle  at  high  velocities 
(2200  fps)  and  temperatures  (16750°  F)  in  which  the  exhaust  jet  turbulent 
length  scale  and  convection  speed  were  measured.  Additionally,  the  laser  was 
successfully  used  to  determine  regions  of  strong  correlation  between  the 
velocity  field  and  the  acoustic  far-field  pressure  for  a  coannular  plug  noz¬ 
zle.  These  recent  results  indicate  that  the  technology  is  now  available  for 
performing  systematic  and  detailed  jet  noise  source  location  diagnostic 
studies  for  jet  nozzles  operating  at  realistic  velocity  and  temperature  condi¬ 
tions  . 
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3.0  GENERAL  REMARKS  ABOUT  THIS  REPORT 


r 


During  the  course  of  the  last  several  years  the  General  Electric  Company 
has  been  engaged  in  a  number  of  theoretical  and  experimental  investigations 
aimed  at  developing  a  better  understanding  and  quantification  of  jet  exhaust 
noise  of  simple  circular  nozzles  and  complex  mechanical  type  exhaust  nozzle 
suppressors.  The  support  for  this  work  has  come  from  a  number  of  Government 
agencies  (FAA/DOT,  Air  Force,  and  NASA)  as  well  as  from  General  Electric 
research  and  development  funds.  Both  scientists  and  engineers  at  the  General 
Electric  Company,  as  well  as  engineers  at  a  number  of  other  industrial  insti¬ 
tutions,  universities,  research  laboratories,  and  Government  research  centers 
have  exerted  and  shared  considerable  time,  talent  and  experience  to  solve  and 
mitigate  the  aeroacoustic  aspects  of  high-velocity  and  temperature  exhaust- 
jet  noise. 

Typically,  one-millionth  of  the  total  jet  flow  power  is  radiated  as 
acoustic  power,  hence  the  difficulty  to  quantify  and  reduce  jet  noise. 
Realistic  exhaust  jets  operate  at  high  velocity  and  temperature.  Therefore, 
complex  and  precise  measurement  and  test  facilities  are  required.  The  jet 
noise  producing  agents  are  aerodynamic  in  nature,  and  are  airborne  -  hence 
rugged,  but  sensitive  real-time  aerodynamic  and  acoustic  measurements  of 
micro  scale  jet  aerodynamic  properties  are  necessary.  Underscoring  these 
problems  is  the  fact  that  theoreticians  are  still  probing  for  the  most  exact 
or  best  theoretical  acoustic  formulations  for  jet  noise  generation  and  for 
acceptable  approximations  for  their  solution.  Considerable  progress  has  been 
made  in  acoustic  testing,  theoretical  acoustic  modeling  and  engineering 
design  procedures  such  as  that  reported  in  the  various  task  reports  under 
this  program  (Report  Series  FAA-RD-76-79) . 

One  perplexing  task  in  verifying  the  "exact"  nature  of  jet  noise  genera¬ 
tion  is  to  develop  the  proper  measurement  tools,  particularly  those  associ¬ 
ated  with  measuring  the  aerodynamic  real-time  flow  properties.  This  report 
discusses  work  carried  out  by  the  General  Electric  Company  to  develop  a  real¬ 
time  velocity  measurement  system,  the  laser  velocimeter. 


Section  4.0  of  this  report  contains  descriptions  and  discussions  of  the 
fundamentals  of  laser  velocimetry,  and  its  unique  problems.  Section  5.0 
contains  an  up-to-date  description  of  General  Electric's  laser  velocimeter 
system,  including  some  unique  real-time  velocity  processing  features  cur¬ 
rently  not  available  in  the  open  literature.  The  theoretical  aeroacoustic 
motivation  for  performing  real-time  velocity  measurements  is  discussed  in 
Section  6.0,  where  a  perspective  is  given  for  performing  real-time,  two- 
point,  space-time  measurements  for  the  purpose  of  measuring  jet  global 
turbulent  structure  properties,  and  for  setting  a  frame  work  for  noise  source 
location  studies.  Sections  7.0  and  8.0  contain  the  theoretical  statistical 
considerations  necessary  for  performing  real-time  velocity  measurements  and 
cross-correlations  with  the  discontinuous  output  of  the  laser  velocimeter. 
Particular  attention  is  devoted  in  these  sections  to  define  the  concepts 
developed  at  GE  which  successfully  permitted  accurate  measurement  of  mean  and 
turbulent  velocities,  turbulent  spectra,  and  two-point  space-time  cross¬ 
correlation  measurements  with  two  lasers,  or  with  one  laser  and  an  acoustic 
microphone.  The  concepts  discussed  in  Sections  7.0  and  8.0  will  provide  a 
basis  for  future  investigators  to  assemble  the  type  of  computational  algorithms 
necessary  for  processing  real-time  data  samples  in  an  economic  fashion. 

Section  9.0  describes  the  laser  velocimeter  test  arrangement  used  in  performing 
many  of  the  reported  advanced  LV  diagnostic  measurements.  Section  10.0 
contains  the  results  of  a  number  of  key  I.V  jet  plume  studies,  including: 
discussions  of  some  high  lights  of  mean  and  turbulent  velocity  measurements 
performed  on  high  velocity  and  temperature  shocked  and  shock  free  conic 
nozzle  flows,  LV  measured  turbulence  spectra,  a  comprehensive  series  of  two- 
point  space-time  measurements  on  a  conic  nozzle  and  an  inverted  flow  coannular 
plug  nozzle  for  turbulent  length  scale  and  convection  velocity,  and  a  series  of 
LV  to  farfield  acoustic  cross-correlation  measurements  on  a  coannular  plug 
nozzle.  The  emphasis  of  all  experiments  was  to  illustrate  the  capability  for 
performing  diagnostic  type  velocity  measurements  in  realistic  high  velocity 
(supersonic)  high  temperature  nozzle  flows.  The  measurements  reported  are 
believed  to  be  unique  and  important  for  future  investigations  of  this  kind. 
Section  11. 0  contains  a  summary  and  recommendations  for  future  work  in  uti¬ 
lizing  the  laser  velocimeter  for  jet  noise  diagnostics. 


4.0  THE  LASER  VELOCIMETER  AND  APPLICATIONS  TO  JET  NOISE  DIAGNOSTICS 


This  section  reviews  the  need  for  and  the  fundamentals  of  laser  veloci- 
metry,  the  historic  background  of  the  technique,  and  unique  problems  and  solu¬ 
tions  in  the  application  of  laser  velocimetry  to  jet  noise  diagnostics. 

4.1  THE  NEED  FOR  A  NON-INTRUSIVE  DIAGNOSTIC  METHOD 

Historically,  in  diagnostics  of  ambient,  subsonic  jets,  a  variety  of 
instrumentation  has  been  available  which  creates  relatively  little  disturbance 
on  the  local  or  global  jet  properties.  Thus,  miniature  pitot-static  probes, 
very  small,  streamlined  microphones,  hot  wire  anemometers,  etc.  were  adequate. 
Such  has  not  been  the  case  for  a  hot,  supersonic  jet.  Various  non-disturbing 
techniques  using  optical  effects  such  as  Schlieren  photography,  holography, 
correlated  two-beam  refraction  and  laser  velocimetry  have  been  attempted  by 
industrial  and  educational  researchers  throughout  the  world.  Ten  years  of 
experience  has  shown  that  the  most  viable  diagnostic  techniques  must  be  accu¬ 
rate  and  highly  localized.  This  accounts  for  the  outstanding  acceptance  of 
laser  velocimetry. 

4.2  HISTORICAL  SKETCH  OF  LASER  VELOCIMETRY 

The  concept  of  measuring  the  velocity  of  a  gas  flow  by  the  Doppler  shift 

of  light  has  been  known  for  a  number  of  years.  In  1948,  Messrs,  Bundy  and 

Strong  of  the  General  Electric  Research  and  Development  Center  published  a 
(8) 

report'  on  the  measurement  of  gas  velocity  in  a  rocket  exhaust,  using  the 

Doppler  shift  of  light  emitted  by  a  trace  of  vapor  of  an  alkali  metal  seeded 

in  the  flame.  The  concept,  however,  was  limited  to  velocity  measurements  of 

hot  gases  where  thermal  excitation  is  sufficient  to  cause  the  alkali  metal 

seeding  material  to  radiate  at  its  resonance  spectral  line.  In  the  1960's, 

the  invention  of  the  laser  provided  a  light  source  with  enough  output  and 

coherence  to  make  heterodyning  techniques  feasible.  In  1964,  Yeh  and  Cummins 

(9) 

also  published  material  on  this  subject  .  However,  improvements  in  utility 
and  arrangement  were  needed.  In  1969,  a  particular  arrangement  of  the  laser- 
powered  Doppler  technique,  called  "Differential  Doppler"  was  first  described 
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by  C.M.  Penney  of  the  General  Electric  Co.  Research  and  Development  Cen- 
ter^1^.  The  geometry  of  this  arrangement  is  such  that  measurements  of  fluid 
flow  velocities  can  be  made  in  situations  where  there  is  access  at  only  one 
side  of  a  stream.  Further,  the  differential  Doppler  arrangement  is  easy  to 
align  and  calibrate.  Today  the  vast  majority  of  users,  worldwide,  are  employ¬ 
ing  the  Differential  Doppler  method,  rather  than  the  original  reference  beam 
arrangement . 

4.3  FUNDAMENTALS  OF  LASER  VELOCIMETRY 

The  laser  velocimeter  is  a  noncontact  air-speed-measuring  instrument  which 
consists  of  a  laser,  seed  particles,  mirrors,  lenses  and  beamsplitters,  a 
photodetector  and  a  signal  processor.  Figure  4-1  is  a  schematic  view  of  the 
basic  differential  Doppler  arrangement.  The  laser  beam  is  split  into  two 
equal  power  beams  that  are  focused  and  crossed  in  space,  where  the  air  speed 
is  to  be  measured.  Interference  "fringes"  are  formed  here  (they  can  be  seen 
with  a  microscope  and  some  smoke).  The  fringes  are  plane,  parallel,  light 
and  dark  sheets  that,  on  edge,  resemble  a  picket  fence.  A  single  dust  particle 
(seed)  passing  through  this  region  scatters  light  in  a  sinusoidal  burst, 
detected  by  a  phototube.  This  is  analogous  to  running  alongside  a  picket 
fence,  rapping  it  with  a  stick.  The  frequency  of  the  sine  burst  is  measured 
by  a  signal  processor.  Figure  4-2  is  a  simplified  schematic  diagram  of  the 
measuring  process.  A  "counter-timer"  is  shown  here,  but  there  are  other 
techniques.  Figure  4-3  is  a  cross  section  to  show  how  the  basic  components 
might  be  arranged  for  a  single  component,  back-scatter  system.  Back-scatter 
means  that  only  a  small  amount  of  the  particle's  scattered  light  is  collected, 
and  in  the  back  direction. 

4.3.1  The  Differential  Doppler  Method 

The  LV  technique  can  take  the  form  of  one  of  two  optical  configurations: 
reference-scatter  (direct  Doppler),  Figure  4-4;  or  differential  Doppler  (dual 
scatter).  Figure  4-5).  The  ease  of  alignment  and  usage  plus  the  inherent 
increase  in  a  signal-to-noise  ratio  have  made  the  differential  Doppler  con¬ 
figuration  predominant  in  the  literature  and  the  method  chosen  by  General 
Electric . 


Figure  4-1.  Schematic  View  of  Laser  Velocimeter 


Figure  4-2.  Schematic  Diagram  of  an  Elementary  LV  Measuring  Process. 
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Figure  4-4.  Reference-Scatter  or  Direct  Doppler  Arrangement. 
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Figure  4-5.  Dual-Scatter  or  Differential  Doppler  Arrangement. 
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An  analysis  of  the  differential  Doppler  method  can  be  understood  by 
considering  the  interference  pattern  created  by  two  intersecting  beams  of 
coherent  radiation,  each  having  a  wavelength  XQ.  Referring  to  the  scattering 
volume  formed  by  the  cross  beams  (Figure  4-6),  constructive  interference 
occurs  when  the  optical  path  length,  A,  between  the  two  beams  differs  by  an 
amount  nX0/2  for  n  odd.  Similarly,  destructive  interference  occurs  when 
A  =  nA0/2  for  n  even.  Planes  of  constant  optical  path  length  difference 
between  the  two  beams  are  parallel  to  the  plane  that  bisects  the  two  beams. 
Thus,  the  fringe  spacing,  S,  is  given  by: 


As  the  particle  passes  through  these  fringes  at  velocity  v,  light  is  scattered 
only  when  it  is  in  a  bright  fringe;  hence,  the  scattering  frequency,  Av,  is 
simply: 


Av 


Aw 

2tt 


=  Vy/S 


or 

Av 


2  Vv  sin  «/2 
lo 


4-2 


where  Vy  is  the  velocity  component  perpendicular  to  the  fringes.  The  inter¬ 
ference  pattern  considerations  are  easier  to  visualize  when  analyzing  the  LV 
as  a  velocity-measuring  device,  since  the  pattern  can  be  thought  of  as  a  light 
grid  placed  in  the  flow  to  define  an  accurate  displacement  increment,  S.  Then, 
tne  local  velocity  normal  to  the  fringes  is  simply  S/t,  where  t  is  the  transit 
time  between  bright  fringes. 


4.3.2  LV  Spatial  Resolution 

The  size  of  the  scattering  volume  is  an  important  consideration  in 
designing  the  LV,  since  it  determines  the  dimension  over  which  any  velocity 
fluctuations  and  gradients  in  the  flow  are  averaged  in  the  measurement.  This 
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Figure  4-6.  Interference  Fringes  Formed  by  Two  Coherent  Crossed  Beams. 


volume  is  determined  by  the  optics  of  the  transmitter  and  the  collector  lenses 
and  by  the  size  of  the  laser  beam.  In  most  LV  arrangements,  the  beams  are 
focused  to  spots  in  order  to  minimize  the  scattering  volume.  Under  these 
conditions,  the  spot  size  is  then  determined  by  the  diffraction  limit  of  the 
focusing  lens. 

Consider  a  beam  of  parallel,  monochromatic  light  of  wavelength  X  with 
uniform  intensity,  passing  through  a  lens  of  focal  length  f.  At  the  focal 
point  of  the  lens,  a  diffraction  pattern  is  formed,  and  the  diameter,  d,  of 
the  central  bright  spot  (called  the  Airy  disc)  is  given  by  the  relation¬ 
ship: 


d  =  1.22  Xf/D  4-3 

where  D  is  the  incident  beam  diameter.  This  same  diffraction  limit  in  the 
size  of  the  spot  at  the  focus  of  the  transmitter  lens  determines  the  scatter¬ 
ing  volume  with  a  laser  beam. 

For  a  laser  beam  incident  on  a  lens  of  focal  length  f>p,  the  intensity  of 
the  laser  beam  is  not  uniform  but  rather  Gaussian  in  character,  and  the  wave- 
fronts  are  normally  not  plane,  but  spherical.  Referring  to  Figure  4-7,  when 


Figure  4-7.  Geometry  of  a  Focused  Laser  Beam. 


a  beam  of  this  nature  (with  diameter  D  at  the  lens  face)  is  incident  on  the 
lens,  the  intensity  in  the  vicinity  of  the  focal  point  is  given  by  (see 
Reference  12): 
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where 

7TX  D2 

U  =  - 2  ’ 

2X  fj 

v  = 

Xfx  ’ 

and  y  and  x  are  coordinates  measured  perpendicular  to  and  along  the  optical 
axis  of  the  beam,  respectively,  with  origin  at  the  focal  point  of  the  lens. 
Thus,  the  focused  beam  is  seen  to  be  a  distorted  ellipsoid  of  revolution  with 
the  intensity  falling  off  sharply  as  a  function  of  x  and  y.  For  the  intensity 
to  drop  to  one  percent  of  its  incident  value,  I0,  the  size  of  the  spot  is: 


=  1.184  Xfx/D  and 

4-5a 

=  7.08  X(fT/D)2. 

4-5b 

These  two  equations  give  the  relative  dimensions  of  the  full  scattering  volume 
without  the  use  of  apertures  and  stops  in  the  collector  optics. 

For  the  differential  Doppler  method,  the  two  beams  intersect  at  the  focal 
point  of  the  transmitter  lens,  and  the  interference  fringes  are  formed  only 
where  the  beam  intersect.  Referring  to  Figure  4-8,  when  the  two  beams  (of 
diameter  D  separated  by  a  distance,  b,  are  focused  by  the  transmitter  lens, 
the  width  of  the  spot  d^  is: 

dD 

cos  (a/2)  ’ 

where  a  is  the  angle  between  the  beams,  and  d0  is  calculated  by  Equation  4-5a. 
Similarly,  neglecting  the  beam  divergence  in  the  vicinity  of  the  focal  point, 
the  length  of  the  intersection,  1°,  is: 
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Figure  4-8.  Transmission  Lens  Scattering  Spot  Size  for  the  Differential  Doppler  Arrangement 
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O'  ~  - 2 - 

°  sin  (a/2) 

As  previously  mentioned,  stops  and  apertures  in  the  collector  optics  are 
normally  used  to  limit  the  size  of  the  scattering  volume  to  a  "data  probe 
volume,"  since,  generally,  the  dimension  of  the  scattering  volume  along  the 
optical  axis,  ZQ  or  Z°,  is  too  long  for  good  spatial  resolution  in  the  flow. 

A  simple  means  for  reducing  this  dimension  by  using  the  collector  optics  is 
shown  in  Figure  4-9,  where  a  narrow  slit  of  height  h  is  placed  in  the  image 
plane  of  the  collector  lens  with  focal  length  fc.  With  this  arrangement,  the 
length  of  the  data  probe  volume,  SI,  in  the  focus  of  the  transmitting  lens  that 
is  seen  by  a  detector  placed  behind  the  slit  is: 

*  =  , 
sin  <p 

where,  p  is  shown  in  Figure  4-9,  and  <j>  is  the  main  scattering  angle.  By  use 
of  the  lens  formula,  one  can  obtain: 

Z  ~  (a/fc  -  1)  h/sin  <p  .  4-8 

Thus,  the  length  of  the  scattering  volume  can  be  controlled  by  a  judicious 
choice  of  the  collector  optics  and  the  scattering  angle. 

With  respect  to  limiting  spatial  resolution,  it  is  possible  to  further 
reduce  the  length  of  the  data  probe  volume  by  discriminating  the  amplitude  of 
the  scattered  light.  Only  when  the  seed  particle  passes  the  middle  of  the 
scattering  volume  is  the  incident  light  (and  therefore  the  scattered  light) 
at  a  maximum.  On  either  side  of  this  point,  the  peak  signature  amplitude 
decreases.  Furthermore,  at  either  end  of  the  intersection  volume  the  signa¬ 
ture  becomes  dual-peaked,  which  can  be  responsible  for  Doppler  dropout  where 
"missing"  electronically  undiscriminated  pulses  will  improperly  add  to  the 
average  Doppler  period.  Thus,  by  amplitude  window  thresholding  (for  particles 
of  uniform  size)  and  use  of  a  dropout  monitor,  the  spatial  resolution  can  be 
increased  beyond  that  generated  by  the  collector  optics. 
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To  summarize,  the  spatial  resolution  of  the  I.V  is  determined  by  both  the 
transmitter  and  the  collector  optics  as  well  as  by  the  use  of  I.V  signal  pro¬ 
cessor  monitors.  The  width  of  the  spot  is  normally  specified  by  the  diffrac¬ 
tion  limit  of  the  transmitter  lens  and  is  given  by  Equation  4-5  or  4-6.  The 
spot  length  is  determined  by  the  collector  optics  containing  appropriate 
apertures  and  stops  to  limit  the  region  over  which  scattered  light  is  collected. 
Dropout  and  amplitude  threshold  windows  also  can  be  provided  as  signal  moni¬ 
tors  to  further  reduce  the  data  volume. 

4 . 4  UNIQUE  REQUIREMENTS,  PROBLEMS  AND  SOI  UTKINS  IN  THE  APPLICATION 

OF  LASER  VELOCIMETRY  TO  JET  NOISE  DIAGNOSTICS 

The  laser  velocimeter  (LV)  is  used  for  the  purpose  of  accurately  sampling 
the  time  history  of  gas  velocity  at  known  points  in  the  nozzle  exhaust  plumes, 
and  in  a  known  direction.  This  sequence  of  velocity  samples  can  be  correlated 
to  near-simultaneous  samples  of  velocity  at  some  other  plume  location,  and 
can  also  be  correlated  with  near  and  far-field  sound  pressure  level  measure¬ 
ments. 

Because  these  diagnostic  LV  data  must  have  sample  time  preserved  as  well 
as  represent  instantaneous  and  localized  velocity  vector  components,  unique 
requirements  are  placed  on  the  LV  system.  Added  constraints  and  problems 
result  when  the  LV  system  is  expected  to  work  in  the  non-ideal  environment  of 
heated  supersonic  jet  exhausts.  The  noise  level  and  temperature  extremes  in 
the  vicinity  of  an  exhaust  nozzle  (model  scale  or  engine  size)  may  seem  incom¬ 
patible  with  a  precision,  laser  powered,  optical  device.  This  section  will 
attempt  to  prove  the  contrary,  that  each  of  the  application  problems,  when 
isolated  and  fully  understood,  can  be  controlled.  Thus,  the  discussion  will 
include:  particles  as  flow  tracers,  use  of  laser  light,  control  volume  size 
and  optical  working  distance,  signal  processing  choices,  angular  acceptance  in 
turbulent  flow,  and  noise  and  other  interfering  effects. 

4.4.1  The  Particle  as  a  Tracer 

The  clean  exhaust  jet  will  not  generally  scatter  the  laser  light  effi¬ 
ciently.  An  external  source  is  usually  required  to  add  tracer  particles 
to  the  gas.  Several  constraints  are  plated  on  these  tracer  particles: 
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1.  must  be  small  in  size  (one  micron  or  less)  to  avoid  "slip"  errors 
in  high  acceleration  fields 

2.  must  be  a  good  light  scatterer 

3.  must  survive  2,000°  R 

4.  must  be  easily  dispersed  in  gas  stream 

4.4.2  Laser  Light  Requirements 

The  light  source  must  be  a  laser  of  high  angular  and  temporal  stability. 
A  complex  set  of  tradeoffs  must  be  made  to  determine  laser  power  required. 

The  following  parameters  or  sub-groups  of  parameters  are  involved: 

1.  laser  wavelength 

2.  laser  coherence 

3.  laser  beam  focusing  optics 

4.  particle  size,  shape,  and  scattering  optical  properties 

5.  scattering  angle  selected 

6.  receiver  collector  solid  angle 

7.  photodetector  quantum  efficiency,  frequency  response  and  noise 
properties 

8.  signal  processor  bandwidth  and  noise  properties 


4.4.3  Control  Volume  Size  and  Optical  Working  Distance 

The  control  volume  generally  is  an  ellipsoidal  shaped  volume,  whose 
dimensions  are  controlled  by  choices  in  the  optics  (see  Section  4.3.2,  par¬ 
ticularly  equations  4-5  and  4-6).  The  optical  working  distance  is  very 
important  in  jet  plume  mapping,  as  it  is  important  to  sample  the  plume  with¬ 
out  exposing  the  LV  head  to  the  jet  blast  itself.  Thus  the  working  distance 
selected  will  depend  on  the  plume  width.  The  control-volume  selected  will  be 
based  on  expected  velocity  gradients  (primarily  that  found  in  the  small  tur¬ 
bulent  eddies  in  the  early  mixing  region  of  the  plume).  A  tradeoff  exists 
here,  because  a  very  small  control  volume  well  suited  to  small  eddy  mapping 
may  lead  to  refraction  problems.  (This  is  discussed  under  noise  and  other 


Various  types  of  LV  signal  processors  may  be  assessed  for  their  measure¬ 
ment  capabilities  of  the  local  fluid  velocity  by  deducing  the  Doppler  frequency 
of  scattered  light  from  a  seed  particle  moving  at  the  fluid  velocity.  These 
signal  processing  designs  include  the  following  major  concepts: 

1.  Spectrum  analyzer 

2.  Doppler  frequency  tracker 

3.  Counter-timer 

4.  Photon  correlator 

5.  Filter  bank 

Early  in  the  history  of  the  LV,  workers  used  a  spectrum  analyzer  for 

processing  of  the  LV  signal.  Such  analyzers,  as  were  and  are  now  available 

in  well  equipped  labs,  employ  a  slow  frequency  scanning  principle.  This  is 

equivalent  to  turning  a  narrow  band  receiver  across  a  frequency  range.  This 

type  of  processor  is  unsuitable  for  measuring  time  dependent  flow  properties. 

An  alternative  approach  to  the  scanning  spectrum  analyzer  is  to  digitally 

record  each  Doppler  burst,  and,  at  some  later  time,  analyze  on  a  digital 

(13) 

computer.  Such  a  technique  was  demonstrated  by  Asher,  Scott  and  Wang 
Advances  in  microcomputers  will  soon  make  on-line  analysis  feasible  and  cost 
effective. 


4. 4. 4.1  Doppler  Frequency  Trackers 

Doppler  frequency  trackers  have  been  in  widespread  use  since  the  late 
1960,s^^,^‘l!t’  .  Commercial  equipment  is  available,  for  example,  from  TSI, 

Inc.  and  DISA,  Inc.^^  Trackers  operate  by  attempting  to  match  the 

frequency  and  phase  of  an  internal,  controllable  oscillator  to  the  incoming 
Doppler  burst.  A  critical  parameter  in  all  trackers  is  the  capture  range. 
Successful  tracking  can  take  place  only  when  the  incoming  Doppler  burst  fre¬ 
quency  is  within  a  certain  ratio  of  the  previous  burst  frequency.  This  ratio 
is  called  the  capture  range.  Usually  this  range  is  0.9  to  1.1.  Low  seeding 
density  and  high  turbulence  in  a  mixing  jet  can  cause  successive  bursts  to 
differ  by  ratios  of  0.3  to  3.0.  Consequently,  the  tracker  is  generally  not 
usable  in  jet  noise  diagnostics. 
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4.4. 4.2  Counter-timer 

Counter-timer  signal  processors  use  conceptually  simple  principles,  but 

have  features  well  suited  to  LV  measurements  in  high  turbulence  gas  flow.  The 
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earliest  presentation  of  the  concept  was  by  W.B.  Jones v  .  Figure  4-2  is  a 
simplified  diagram  of  Jones'  arrangement.  The  high  frequency,  brief  burst 
from  the  photomultiplier  tube  is  high  pass  filtered  to  remove  the  gaussion, 
low  frequency  envelope  (sometimes  referred  to  as  a  "pedestal").  The  signal 
is  then  amplified  and  clipped  and  enters  the  preset  counter  which  outputs  a 
single  pulse  of  duration  Nt,  where  t  is  the  burst  period  and  N  is  the  preset 
number  of  cycles  to  be  counted.  The  pulse  enters  a  width  to  height  converter 
that  has  output  Kt  as  an  amplitude.  Thus  Kt  is  proportioned  to  the  transit 
time  of  the  seed  particle  through  the  N  fringes.  The  transit  time  signal, 

Kt,  may  be  recorded,  or  applied  to  a  display  oscilloscope,  or  both. 

The  counter-timer  processor  is  widely  used  in  low  speed/density  flow 
where  high  precision  is  needed  in  the  measurement.  Advanced  developments 
relating  to  jet  noise  diagnostics  are  described  in  Section  3.0. 

4. 4. 4. 3  Photon  Correlation 

Photon  correlation  methods  and  their  application  to  laser  light  scatter¬ 
ing  techniques  have  been  discussed  by  Pike^1^  and  by  Kalb^^.  The  basis  of 
this  approach  is  the  detection  of  individual  quanta  of  light  energy  (photons) 
recorded  in  fixed  time  intervals.  Using  appropriate  digital  processing,  the 
recorded  photon  "counts"  can  be  shown  to  represent  the  autocorrelation  func¬ 
tion.  This  signal  then  contains  time-averaged  velocity  and  turbulence  inten¬ 
sity  information.  For  this  reason,  a  photon  signal  processor  could  only  be 
used  in  a  very  inefficient  way  in  jet  plume  diagnostic  velocity  measurements. 
It  has  been  calculated  that  a  velocity  spectral  analysis  using  photon  correla¬ 
tion  would  require  data  collection  times  in  excess  of  that  needed  by  the 
counter-timer,  for  example,  by  a  factor  in  excess  of  100. 

4. 4. 4. 4  Filter  Bank 

Filter  bank  signal  processing  for  the  LV  was  first  reported  in  1968  by 
(14) 

Rolfe  et  al.  Although  showing  feasibility  of  using  a  set  of  overlapping 
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bandpass  filters  for  an  efficient  spectral  analysis,  the  arrangement  was 
limited  in  useiulness.  The  arrangement  averaged  many  Doppler  bursts,  finally 
yielding  mean  velocity  and  turbulent  intensity.  Times  of  arrival  of  individual 
bursts  were  not  preserved.  The  potential  for  individual  burst  processing, 
however,  was  later  pursued  at  General  Electric,  and  a  filter  bank  processor 
that  is  suitable  for  jet  plume  time  dependent  measurements  was  developed  and 
applied  as  described  in  Section  5.0. 

4.4.4. 5  Summary 

The  counter-timer  and  the  filter  bank  were  selected  by  GE  to  process  a 
velocity  sample  from  just  one  particle  traverse  of  the  LV  control  volume. 

The  spectrum  analyzer,  Doppler  frequency  tracker,  and  photon  correlator  were 
found  to  be  either  not  capable  of  or  inefficient  in  handling  single  particle 
Doppler  burst  data. 

4.4.5  Angular  Acceptance  in  Turbulent  Flow 

As  the  flow  turbulence  increases,  the  velocity  components  show  increas¬ 
ing  fluctuations.  The  signal  processor  must  contend  with  the  wide  frequency 
range  inherent  in  this  process.  In  addition,  turbulence  brings  about  angular 
trajectory  fluctuations  about  the  jet  centerline  axis.  This  illustrated  in 
Figure  4-10,  tracer  particle  trajectories  in  the  control  volume.  Three 
illustrative  particle  trajectories  of  -60°,  0°,  and  +60°  are  shown,  where  0° 
represents  the  jet  axis.  Sixteen  (16)  fringes  were  placed  in  the  control 
volume  in  this  example.  The  0°  particle  crosses  all  16  fringes,  but  only  8 
fringes  are  crossed  by  the  -60°  and  +60°  particles.  The  governing  equation 
of  the  number  of  fringes  crossed  or.  the  prcbe  control  volume  center  is  a 
cosine  law,  and  is  illustrated  in  Figure  4-11.  As  an  example  of  the  signal 
processing  limits,  suppose  a  control  volume  of  only  ten  (10)  fringes  was 
used,  and  the  signal  processor  required  eight  (8)  good  cycles  (fringes)  of 
burst  signal.  In  Figure  4-11,  the  dashed  line  marked  "A"  would  represent 
the  angular  limit  line  for  this  case,  where  only  0.8  of  the  total  fringes  are 
crossed.  It  may  be  seen  that  line  A  intersects  the  cosine  curve  at  approxi¬ 
mately  37°.  Thus,  particles  beyond  ±37°  would  not  be  measurable.  If  it  was 
determined  that  the  angular  acceptance  of  ±37°  was  not  sufficient,  then  the 
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ratio  of  processor  cycles  required  over  total  fringes  in  the  control 
volume  would  have  to  be  decreased.  For  example,  line  "B"  in  Figure  4-11 
shows  a  ±60°  angular  acceptance  when  1/2  of  the  fringes  are  crossed. 

The  examples  given  are  simplified  somewhat,  in  that  some  particle  tra¬ 
verses  through  the  control  volume  will  be  off-center  and  cross  fewer  fringes 
than  illustrated  above. 

4.4.6  Noise  and  Other  Interf erring  Effects 

This  section  outlines  the  problems  of :  optical  and  electrical  noise, 
interference  from  background  light,  smoke  interference,  optical  refraction, 
and  environmental  protection. 

4.4.6. 1  Optical  and  Electrical  Noise 

Optical  and  Electrical  Noise  stems  from  several  sources:  A  small  par¬ 
ticle  of  non-spherical  shape  may  rotate  as  it  passes  through  the  LV  control 
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volume,  creating  intensity  modulation  in  the  scattered  light  collected 

The  focused  and  crossed  laser  beams  do  not  have  a  perfect  sinusoidal  light 

distribution,  but  have  what  may  be  thought  of  as  secondary  interference  pat- 

(12) 

terns  .  The  photosensor  exhibits  "shot  noise"  that  stems  from  the  very 

nature  of  light:  light  travels  in  packets,  called  photons,  that  vary  in  a 

statistically  predictable  and  sometimes  annoying  manner  -  even  from  a  so-called 
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steady  light  source  .  The  impact  of  these  noise  sources  is,  with  a  single 
particle  in  the  control  volume,  a  range  of  10  to  30  dB  signal  to  noise  ratio 
on  the  Doppler  burst  itself. 

The  most  influential  parameter  is  the  scattered  light  power  incident  on 
the  photo-detector.  Any  method  used  to  increase  the  detected  light  power 
will  improve  the  signal  to  noise  ratio. 

4. 4. 6. 2  Interference  from  Background  Light 

Light  received  from  any  source  other  than  the  scattered  laser  light  from 
the  particle  constitutes  an  interference  source.  Therefore,  hot  particle 
incandescence,  illuminated  background  from  roomlight,  sunlight,  and  even  the 
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laser  beams,  may  degrade  the  signal.  The  laser  line  filter  will  remove  over 
99%  of  unwanted  broadband  sources.  The  filter  will  not  remove  background 
laser  light.  This  light  may  only  be  excluded  by  keeping  surfaces  out  of  the 
control  volume  during  data  acquisition. 

4.4. 6. 3  Smoke  Interference 

Smoke  interference  is  related  to  the  background  light  problem.  The 
laser  light  scattered  from  a  heavy  concentration  of  smoke  particles  presents 
the  same  interference  as  a  surface  in  the  control  volume.  In  addition,  smoke 
between  the  control  volume  and  the  LV  head  may  attenuate  the  signal  by  absorp¬ 
tion  and  scattering. 

4 . 4 . 6 . 4  Optical  Refraction 

Optical  refraction  is  the  bending  of  light  rays  in  the  air  space  between 
the  LV  head  and  the  jet,  and  within  the  jet.  Temperature  gradients  outside 
and  within  the  jet  plume  and  static  pressure  variations  within  the  plume 
create  the  optical  index  gradients  which,  in  turn,  cause  the  refraction. 

Under  severe  refraction,  the  laser  beams  uncross.  The  loss  of  the  interference 
fringes  during  the  uncrossed  condition  stops  the  data.  Partially  uncrossed 
beams  result  in  shortened  Doppler  bursts  and/or  reduced  signal  strength. 

The  immersion  depth  of  the  optical  path  in  the  hot  jet  affects  the  amount  of 
refraction. 

4.4. 6. 5  Environmental  Protection 

Environmental  protection  of  the  LV  head  involves  several  areas.  The 
mechanical  mountings  of  the  laser,  mirrors,  lenses  and  receiver  optics  must 
be  designed  to  withstand  temperature  changes  without  undue  strain.  Use  of 
alloys  of  equal  thermal  expansion  coefficient,  and  high  thermal  conductivity 
reduce  the  strain  problem.  Acoustic  noise  and  vibration  effects  can  be  con¬ 
trolled  by  designing  for  the  highest  possible  resonant  frequency  in  the  LV 
head  and  internal  optics  mountings.  Smoke  and  dust  must  be  prevented  from 
depositing  on  the  optics.  Near-air-tight  construction,  to  reduce  "breath¬ 
ing",  helps.  Positive  purge  with  a  filtered  air  source  is  even  better. 
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5.0  DESCRIPTION  OF  GENERAL  ELECTRIC'S  LASER  VELOCIMETER 


This  section  presents  an  up-date  on  the  development  and  status  of  General 
Electric's  current  arrangement  of  the  real-time  diagnostic  laser  velocimeter 
systems.  The  topics  covered  are: 

•  Laser  velocimeter  optical  head 

•  Signal  processors 

•  Seeding  arrangements 

•  Data  acquisition  and  reduction 

5.1  LASER  VELOCIMETER  OPTICAL  HEAD 

Since  1970,  six  LV  heads  have  been  designed,  built  and  used.  Four  are  cur¬ 
rently  in  use.  Table  5-1  summarizes  the  characteristics  of  these  LV  heads. 

All  are  of  the  single  component,  differential  Doppler  arrangement  (defined  in 
Section  4.0).  All  except  the  Model  10-1  are  of  the  backscatter  arrangement. 

The  optical  working  distance  varies  from  10  inches  to  85  inches.  The  shorter 
distances  are  suitable  for  probing  turbomachinery  blading.  The  longer  optical 
working  distances  have  been  designed  for  probing  jet  exhaust  plumes,  where  no 
physical  part  of  the  head  is  allowed  to  contact  the  high  speed  exhaust  jet. 

In  Table  5-1,  the  column  "Coaxial  or  Offset"  refers  to  whether  the  laser 
beam's  bisector  has  been  arranged  to  be  centered  inside  of  the  received 
cone  of  scattered  light  ("coaxial”)  or  outside  of  this  cone  ("offset") . 

Note  also  that  view  cone  requirements  vary,  and  probe  volume  dimensions  are 
quite  different. 

The  two  LV  heads  used  in  the  cross-correlation  measurements  (discussed 
in  Section  10.0)  are  the  models  85-1  and  85-2,  lines  4  and  6  on  Table  5-1. 
Although  outside  case  dimensions  are  slightly  different,  the  probe  volumes 
were  matched  as  closely  as  possible.  In  two-point  velocity  correlations, 
one  LV  is  operated  on  the  514.5  nm  argon-ion  laser  "green"  line.  The  other 
head  is  operated  on  the  488  nm  "blue"  line.  Each  photomultiplier  tube  is 
filtered  with  the  appropriate  laser  line  filter.  Because  of  the  two  colors 
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Table  5-1.  Laser  Velocimeters  Designed,  Built,  and  Used  at  General  Electric  Company' 
Aircraft  Engine  Group  Measurements  and  Sensor  Development  Lab. 


used,  probe  volumes  can  be  moved  in  close  proximity  to  one  another,  even 
overlapping,  without  interference.  The  general  arrangement  of  the  85-1  and 
85-2  LV  heads  are  shown  in  Figure  4-3.  The  argon-ion  laser  supports  a 
3X  beam  expanding  telescope.  One  beamsplitter  and  three  mirrors  are  used 
to  form  the  two  equal  power  laser  beams  of  about  0.7  watts  each.  The  beam 
expander  is  adjusted  so  as  to  focus  the  laser  beams  about  90  inches  away, 
to  a  spot  diameter  of  0.018  inches.  No  other  lenses  or  optics  are  used  in 
the  transmitter  side  of  the  LV.  On  the  receiver  side,  a  pair  of  symmetrical, 
6-in.  diameter,  acromatic  lenses  are  used  to  collect  a  4°  cone  of  scattered 
light  and  focus  through  a  field  stop  to  the  filter  and  detector.  The  detector 
used  is  a  photomultiplier  tube,  0.75  inch  diameter  "S-20"  cathode  and  potted 
dynode  resistors.  A  wideband  preamplifier  is  located  about  10  inches  from  the 
detector.  Frequency  response  of  the  preamplifier  is  1  kHz  to  150  MHz  +  1  db. 
Voltage  gain  is  100X.  Peak  output  is  ±  1  volt  into  a  50  0  transmission  line 
(such  as  type  RG8  cable) . 

The  LV  box  construction  is  1/2  inch  thick  aircraft  aluminum  alloy  sheet, 
with  a  1/16  inch  thick  plumber's  lead  sheet  bonded  to  the  inside.  The  alumi¬ 
num  alloy  provides  a  high  mechanical  resonant  frequency  and  high  thermal 
conductivity.  The  lead  sheet  provides  a  high  mechanical  damping.  For  opera¬ 
tion  in  proximity  to  scale  model  exhaust,  no  additional  acoustic  protection 
was  found  to  be  necessary.* 

5.2  SIGNAL  PROCESSORS 

Since  1970,  six  signal  processors  have  been  designed,  built  and  used, 
encompassing  five  designs.  Four  are  currently  in  use.  Table  5-2  summarizes 
the  characteristics  of  these  processors.  The  first  four  lines  of  the  table 
are  models  based  on  the  counter-timer  technique.  The  fifth  line  describes 
the  Mark  III  signal  processor,  based  on  the  filter  tank  technique.  A  com¬ 
parative  description  of  these  techniques  is  given  in  Section  4.0. 

*  The  LV  has  been  used  in  diagnostic  testing  of  a  J79  engine  at  full  and 
partial  power.  During  such  testing  a  time  delay  relay  tube  in  the  laser 
powersupply  box  (about  two  feet  from  the  LV  head)  was  broken.  Wrapping  a 
lead-vinyl  acoustic  blanket  around  the  laser  supply  box  solved  this 
problem  and  no  further  damage  to  the  LV  was  experienced  during  engine  testing. 
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To  the  counter-timer  basic  arrangement  (line  1,  Table  5-2)  Mossey 
added  a  wide  frequency  range  pedestal  remover  (line  2),  a  precision  counter¬ 
timer-comparator  (line  3)  and  a  period-to-velocity  converter  (see  Chapter  5 
of  Reference  23,  and  Reference  24).  This  formed  the  basis  of  the  Mark  I  signal 
processor.  Additional  refinements  were  added  by  Asher  (25).  The  success 
of  the  Mark  I  processor  paved  the  way  to  development  of  the  more  elaborate 
Mark  II  processor  (line  4) .  The  Mark  II  processor  has  an  advanced  check 
system.  Two  of  this  type  of  processor  are  currently  in  use.  The  significant 
advantage  of  the  Mark  II  is  the  greatly  reduced  error  rate  in  the  data  out¬ 
put  stream  over  any  previous  processor. 

Because  of  these  discretionary  features  built  into  the  Mark  II  signal 
processor,  the  data  output  rate  drops  as  the  signal  to  noise  ratio  (S/N)  falls 
below  about  20  db.  In  order  to  improve  the  data  output  rate  under  specific 
low  S/N  conditions.  General  Electric  developed  the  Mark  III,  a  completely 
different  type  of  signal  processor,  shown  on  line  5  of  Table  5-2. 

Both  the  Mark  II  and  Mark  III  signal  processors  have  been  used  in  diag¬ 
nostic  and  cross-correlation  velocity  measurements.  A  detailed  operation 
description  of  each  processor  is  given  below. 

5.2.1  Mark  II  LV  Signal  Processor 

The  Mark  II  LV  signal  processor  is  illustrated  in  Figure  5-1.  The 
amplified  Doppler  burst  signal  is  applied  to  the  wideband  pedestal  remover. 

The  square-wave  like  output  is  illustrated  in  Figure  5-1.  This  output  is 
applied  to  a  four  stage  binary  counter,  and  also  to  an  under-frequency  tester. 
The  binary  counter  feeds  eleven  (11)  decoders,  and  each  decoder  has  a  unique 
logical  function.  Each  decoder  feeds  a  time-to-amplitude  converter.  The 
first  nine  (9)  decoders  are  for  specific  single  cycles  within  the  burst. 

All  of  these  decoders,  through  their  respective  time-to-amplitude  converters, 
feed  the  all-cycle  comparator  test.  In  this  test,  one  of  the  cycle  periods 
is  chosen  as  the  reference  period.  Then,  all  other  periods  are  compared 
with  this  reference.  If  any  periods  deviate  from  the  reference  by  more  than 
the  preset  fractional  value,  an  "invalidate"  signal  results.  Assuming  all 
tests  were  "OK",  a  velocity  signal  output  and  a  Data  "OK"  output  result.  The 
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Figure  5-1.  Mark  II  LV  Signal  Processor. 
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velocity  signal  may  be  selected  to  be  calculated  from  four  (4)  cycles  of  the 
burst  or  from  eight  (8)  cycles  of  the  burst,  by  a  manual  switch  setting. 

For  the  LV  diagnostic  tests  performed  to  date,  the  nominal  setting  used 
for  data  validation  is  the  eight  (8)  cycle  setting.  In  addition  to  the 
under-frequency  and  all-cycle  tests,  a  maximum  and  minimum  signal  test  is 
done  on  the  amplitude  of  the  burst.  This  helps  to  qualify  the  particle 
size,  since  scattered  light  is  dependent  to  about  the  fourth  power  of  scat- 
terer  size. 

This  is  believed  to  be  the  most  complete  comparative  validation  tech¬ 
nique  in  existance  on  a  counter-time  processor.  False  data  have  been 
practically  eliminated. 

5.2.2  Mark  III  LV  Signal  Processor 

The  Mark  III  LV  signal  processor  is  illustrated  in  Figure  5-2.  This  pro¬ 
cessor  was  developed  to  produce  the  highest  possible  data  rate  under  poor 
signal  conditions.  The  model  for  the  design  is  the  "optimal  receiver"  -  a 
technique  for  detecting  a  signal  and  estimating  its  frequency  and  arrival  time. 
The  basis  of  the  receiver  is  a  set  of  filters,  spanning  the  frequency  range 
of  Interest,  but  matched  to  the  expected  signal  envelope  shape.  The  analysis 
is  given  by  Whalen  in  Reference  26.  Although  the  principle  has  been  well 
known  in  communications  engineering,  the  application  has  been  largely  avoided 
due  to  its  high  cost  and  complexity.  As  an  example,  to  obtain  a  frequency 
resolution  of  1/2%  by  a  conventional  filter  bank,  the  filters  would  be  spaced 
in  increments  of  1%.  For  a  processor  range  of  0.5  to  65  MHz,  a  total  of  490 
filters  would  be  needed. 

In  1972,  a  technique  was  conceived  that  would  greatly  reduce  the  number 
of  filters  needed  for  a  given  resolution.  In  this  technique,  instead  of 
choosing  the  largest  output  filter,  the  outputs  of  all  filters  with  activity 
over  a  threshold  are  used  in  a  calculation  program  for  an  arithmetic  mean 
frequency.  The  Mark  III  processor  uses  this  innovation,  and  has  a  resolution 
of  1/2%  of  frequency  with  only  120  filters  to  cover  the  0.5  to  65  MHz  range. 

As  in  the  Mark  II,  the  Mark  III  processor  gives  a  complete  data  output 
(velocity  and  arrival  time)  for  each  single  particle  traverse.  The  Mark  III 
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Matched  Filters 
Spanning  Frequency  Range 


Figure  5-2.  Mark  III  LV  Signal  Processor. 
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also  has  tests  that  are  applied  to  identify  "good"  data.  In  Figure  5-2  a 
signal  is  examined  for  level  in  a  "max-min  test".  A  "bandwidth  test"  is  done 
to  establish  the  spectral  distribution  of  the  signal.  For  example  a  very 
noisy  burst  will  result  in  a  wide  bandwidth  or  excess  number  of  filters 
having  activity.  Only  if  all  tests  are  met,  will  a  DATA  "OK"  pulse  result. 

5.3  SEEDING  ARRANGEMENTS 

For  all  the  scale-model  tests  performed,  an  aluminum  oxide  powder  of 
about  0.7  microns  mean  diameter  was  used.  Depending  on  the  model  tested, 
seeding  is  provided  to  the  air  entering  the  combustor  of  the  inner  stream 
flow,  the  outer  stream  air  (for  coannular  flow  nozzles),  and  the  entrained  aii 
Seed  generation  is  provided  by  a  fluidized  bed  approach  which  Vias  a  remotely 
variable  rate  control.  Figure  5-3  is  a  diagram  of  the  generator,  and  Figure 
5-4  is  a  photograph  of  the  seed  generator  and  feed  and  control  lines  as  in¬ 
stalled  at  the  General  Electric  Anechoic  Jet  Noise  Test  Facility.  The 
fluidized  bed  is  an  ideal  device  to  entrain  the  particles  in  the  air.  The 
fluidized  bed  produces  an  updraft  that  produces  a  slow  "boiling"  action  in 
the  bed  which  allows  only  the  finest  particles  to  be  entrained.  Large 
cluster  and  agglomerates  tend  to  stay  behind  with  this  approach.  As  a  re¬ 
sult,  the  entrained  seed  may  be  sent  through  ordinary  tubing  and  pipe, 
including  90°  elbows,  tees,  and  valves  over  long  distances  with  no  clogging 
or  drop-out  problems.  Generally,  a  commercial  free-flow  agent,  such  as  fumed 
silica,  is  added  to  assist  feeding  and  prevent  agglomeration. 


5.4  DATA  ACQUISITION  AND  REDUCTION 

The  full  sensor  system  at  General  Electric's  Anechoic  Jet  Noise  Test 
Facility  consists  of  two  laser  velocimeters,  a  network  of  microphones  (both 
near-field  and  far-field),  two  closed  circuit  television  systems  (one  moni¬ 
tors  the  LV  beams,  the  other  monitors  the  nozzle  under  test),  and  a  complete 
gas  and  metal  temperature,  and  gas  pressure  sensor  network  for  scale  model 
nozzle  operation.  Two  minicomputers  are  used:  one  for  nozzle  operation, 
the  other  for  jet  plume  diagnostic  data  acquisition  and  reduction.  This 
section  describes  only  the  jet  plume  diagnostic  data,  control,  acquisition 
and  reduction.  This  equipment  is  shown  in  schematic  form  in  Figure  5-5. 
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Figure  5-5.  Data  Acquisition  and  Reduction. 
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Individual  controls  are  used  for  each  LV.  Each  LV  has  its  own  original 
processor  assigned.  Each  processor  has  its  own  pulse  height  analyzer  and 
plotter. 


The  General  Electric  LV  processor  outputs  are  also  directly  linked  to  a 
PDP  11/45  minicomputer  system  which  is  used  for  the  following  data  reduction 
and  cataloging  tasks: 

•  Log  nozzle  conditions  and  plume  sensor  locations. 

•  Count  and  log  number  of  LV  data  samples  per  location. 

•  Plot  histogram  of  velocity  samples. 

•  Plot  cross-correlation  coefficients  (in-jet  in-jet  and  in-jet  far- 
field  microphone). 

•  Calculate  and  log  correlation  coefficients,  mean  and  root-mean- 
square  values,  confidence  intervals,  histogram  data. 

The  minicomputer  lias  a  multiple  input  analog  to  digital  converter, 
floating  point  processor,  core,  disc  and  bulk  magnetic  tape  storage,  keyboard 
control,  storage  scope  display  and  hard  copy  unit.  The  data  reductions  neces¬ 
sary  to  control  the  experiment  and  confirm  that  adequate  data  has  been 
acquired  (confidence  interval)  arc  performed  within  seconds  of  each  data  col¬ 
lection  interval.  This  allows  the  engineer  to  make  a  decision  for  or  against 
additional  LV  samples  before  actuating  the  LV  system  to  the  next  plume  loca- 
t  ion . 

All  of  the  above  data  is  provided  on-line  during  test,  and  is  all  filed 
in  disk  storage.  The  stored  data  is  used  later  for  data  reduction  for  length 
and  scale  of  turbulence,  convm  tion  velocity,  spectrum  of  turbulence  and 
derivatives  of  cross-correlation  functions.  The  data  reduction  programs  used 
for  on-line  data  analysis  and  reduction  of  the  jet's  turbulence  properties 
are  discussed  in  Section  7.0  and  8.0. 


6.0  THEORETICAL  AERO-ACOUSTIC  FRAME  WORK  FOR  REAL  TIME  MEASUREMENTS 


The  problem  of  performing  meaningful  turbulent  flow  measurements  for  the 
purpose  of  better  understanding  jet  noise  continues  to  be  a  difficult  one  for 
technology  and  physics  in  general.  Until  recently,  the  prospect  of  performing 
real  time  velocity  measurements  in  realistic  jet  flows  was,  practically  speak¬ 
ing,  impossible.  With  the  advent  of  today's  laser  velocimeter,  the  prospect 
has  changed  and  moved  in  a  direction  of  high  feasibility.  It  is  only  now 
that  the  acoustic  design  engineer  has  at  his  disposal  the  type  of  advanced 
instrumentation  necessary  for  performing  the  type  of  detailed  measurements 
needed  to  better  understand  and  postulate  more  accurate  aero-acoustic  models. 

Jet  noise  generation  and  reduction  is  directly  linked  to  the  aerodynamic 
nature  of  turbulent  jets  in  the  form  of  what  is  referred  to  as  the  "aerodynamic 
source  team".  The  definition  of  this  term  can  be  analytically  described  in 
terms  of  direct  correlations  between  the  turbulent  velocity  and  the  acoustic 
far-field  pressure,  or  in  terms  of  turbulent  structure  properties  such  as 
turbulent  length  scale,  urbulent  convection  speed,  and  turbulence  spectra  of 
the  moving  eddies.  The  latter  turbulence  properties  can  be  measured  by  so- 
called  turbulent  velocity  two-point  cross  correlation  measurements.  Some  of 
the  governing  equations  which  are  used  for  describing  the  aerodynamic  nature 
of  jet  noise  are  discussed  below. 

6.1  GOVERNING  EQUATIONS 

The  generation  of  aerodynamic  noise  by  free  turbulence  received  its  first 

(27  28 ) 

quantitative  formulation  by  Lighthill.  ’  1  The  principle  quantitative  step 
by  Lighthill  was  his  rearrangement  of  the  continuity  and  momentum  equations 
to  yield  a  wave  equation. 

For  the  continuity  equation  we  have: 

^  +  .  o  (6-1) 

■)t  dXf 
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where  the  coordinate  system  is  defined  in  the  sketch  on  the  following  page: 


The  subject  of  applying  the  approximations  and  extensions  to  Lighthill's 

original  work  are  numerous.  Some  of  GE's  recent  publications  are  contained 

in  References  7,  23,  29,  30,  31,  32,  and  33.  However,  a  team  of  researchers 

at  the  University  of  Toronto  under  H.S.  Ribner  has  been  performing  analysis 

( 3  A  33 

and  measurements  to  obtain  a  direct  correlation  of  noise  and  flow.  ’  ’ 

36  37) 

’  From  the  point-of-view  of  a  direct  theory  experiment  verification, 

and  with  an  objective  of  learning  how  to  use  direct  noise  source  location 
measurements  to  help  in  the  design  of  jet  noise  suppressors,  the  work  of  the 
Toronto  team  is  unique  and  attractive.  A  summary  of  useful  concepts  for 
applying  in-jet  type  measurements  to  jet  noise  diagnostics  is  presented 
below. 


6.2  TYPES  OF  AERO-ACOUSTIC  MODELS  FOR  JET  NOISE  DIAGNOSTICS 
6.2.1  Indirect  Relationships 

To  obtain  estimates  for  acoustic  intensity,  a  considerable  simplifica¬ 
tion  to  solving  Equation  6-4  results  in  expressing  the  Lighthill  stress 

(38) 

tensor  contraction  in  the  Proudman'  form. 


x.x. 
i  1 
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(37) 

W.T.  Chu  has  shown  that  the  spectrum  of  noise  per  unit  volume  of 
turbulence  may  be  estimated  using  the  Proudman  approximation: 


1  (x.  f) 


T  a  T 
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u  5/2  4 
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(6-5) 


where 


L  , 
x  f 


Typical  length  scale  and  frequency  of  turbulence 
Radian  frequency  in  the  turbulence 


0  = 


oj2  l2 

(1-M  Cos  0)2  +  (-V21) 


1/2 


Convection  Mach  Number  of  the  Turbulent  eddy 

(Vc/a0) 


T  =  f’o  Uv 

XX  Ox 


Ux  =  The  instantaneous  velocity  component  in  the  x-direction 

U  +  u'  (t) 
x  x 


Thus,  from  Equation  6-5,  note  that,  if  the  point  properties  of  the  turbulence 

(L  ,  U',  U  ,  Wv )  can  be  measured,  an  estimate  of  the  far-field  acoustic  radi- 
x’  x  x  L 

ation  per  unit  volume  of  jet  can  be  made.  Sections  7  and  8  below  describe 
the  algorithms  necessary  to  obtain  these  real-time  measurements  using  the 
laser  velocimeter.  Suffice  to  say  here  that  the  mean  velocity  and  turbulent 
velocity  can  be  determined  by  a  single  point  measurement.  However,  to  deter¬ 
mine  the  convection  speed,  length  scale  of  turbulence,  and  the  turbulence 
spectra  of  the  convecting  eddy  requires  the  use  of  two-point  space-time  cross 
correlation  measurement  technology. 


r.4 


6.2.2  Direct  Relationships 

In  formulating  more  direct  relationships,  or  correlation  between  the  jet 

flow  and  noise,  a  number  of  investigators  have  formulated  unique  techniques 

which  can  be  used  for  various  applications  (see  References  35,  36,  37,  39, 

40,  and  41).  Some  of  the  reported  work  was  formulated  for  correlations  of 

two  pressure  signals  (References  39,  40,  and  41)  and  others  are  formulated  as 

a  direct  cross-correlation  between  the  source  velocity  field  and  the  acoustic 

radiation  field.  In  particular,  the  work  of  Lee  and  Ribner  ,  and  Seiner 
(42) 

and  Reethof  are  notable.  Basically,  however,  all  the  techniques  advanced 
were  limited  in  their  demonstration  to  subsonic  ambient  jets  due  to  instru¬ 
mentation  limitations.  The  laser  velocimeter  does  not  have  the  limitations 
of  these  older  measurement  techniques  and  thus  the  LV  was  adapted  and  de¬ 
veloped  for  the  direct  causality  type  measurements. 

To  obtain  a  working  relationship  relating  the  turbulent  jet  velocity 
field  to  the  far-field  acoustic  pressure  field,  use  of  Proudman's  form  of 

(35) 

Lighthill's  integral  for  aerodynamic  noise  is  used  again  .  An  illustrative 
equation  is  an  expression  of  acoustic  intensity  at  x  due  to  a  unit  volume  of 
turbulence  at  y. 
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where 


'  2  i 

ux  P 


=  u2  (y)  p'  (X;  y,  t ) 


=  Retarded  Time  (x/a0) 


=  U  -  u'  (t) 

X  X 


Hence,  once  the  value  of  R  2  ,  and  the  retarded  time  t-  are  known,  in  prin- 

ux  P  ° 

ciple  the  acoustic  intensity  of  the  noise  can  be  determined.  Although  the 
concept  is  simple  in  principle,  carrying  out  the  proper  measurement  using  a 
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laser  velocimeter  is  quite  complicated.  The  theoretical  considerations  and 
computational  algorithms  used  in  General  Electric's  I.V  system  are  described  in 
Sections  7.0  and  8.0. 

One  additional  concept  which  has  been  introduced  into  6-6  is  the  decom¬ 
position  into  two  types  of  turbulence  noise  sources.  To  do  this,  the  instan¬ 
taneous  turbulent  velocity  field  Ux(y,t)  can  be  written  as: 

Ux(y,t)  =  Ux  +  ux(y,t) 

The  cross-correlation  function  Ru2pi  can  be  expressed  in  terms  of  two  cross¬ 
correlations  as  follows: 


u2p'  =  E  ju2p'j 

=  E  j(02  +  2  Ux  u’x  +  u'2)  p-| 

=  u2e(p’)°  +  2UxE(u'xp')  +  E  (u'2  p') 


where  E  {•}  denotes  statistical  expectation. 

The  first  term  on  the  right-hand  side  of  Equation  6-7  may  be  identified  as 
the  shear  noise  cross-correlation  term,  while  the  second  term  may  be  iden¬ 
tified  as  a  self  noise  cross-correlation  term.  Later,  the  magnitude  of  these 
measured  terms  will  be  performed  and  some  interesting  conclusions  will  be 
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drawn. 


7.0  THEORETICAL  CONSIDERATIONS  FOR  REAL  TIME  MEASUREMENTS 
USING  THE  LASER  VELOCIMETER. 

The  previous  sections  have  established  the  importance  of  mean  and  turbu¬ 
lence  velocity  measurements,  velocity  histograms,  and  time  series  statistics 
such  as  velocity  spectrum,  and  cross  correlation  functions  to  the  aero-acoustic 
engineer.  This  section  develops  the  theory  necessary  for  making  these  measure¬ 
ments  with  the  laser  velocimeter.  The  reason  that  special  measurement  tech¬ 
niques  must  be  developed  for  the  laser  velocimeter  is  that  the  LV  makes  dis¬ 
crete  measurements  randomly  in  time  (measurements  are  made  only  when  particles 
arrive  at  the  probe  volume).  The  traditional  theories  of  measurement  apply 
only  for  continuous  signals  or  for  signals  sampled  in  equal  spaced  time 
intervals.  New  theories  of  measurement  are  necessary  to  treat  the  random  tim< 
sampling  that  occurs  in  laser  velocimeter  measurements. 

Two  types  of  measurements  (histograms  and  correlation  functions)  can  be 
made  from  which  mean  velocity,  turbulence  velocity,  and  velocity  spectrum 
measurements  can  be  obtained. 


7 . 1  HISTOGRAM  MEASUREMENTS  FROM  LASER  VELOCIMETER  DATA 

A  histogram  is  an  estimate  of  the  first  order  probability  density  function 
of  the  amplitude  of  a  signal.  To  obtain  a  histogram,  the  range  of  the  input 
signal,  x(t),  is  divided  into  classes  bounded  by  values  .  A  counter,  k^,  is 
incremented  for  each  input  if  <x(t)  <x.+^ .  This  is  done  for  N  independent 
samples  of  x(t),  such  that 


where  P^  is  an  estimate  of  the  probability  that  a  sample  of  x(t)  is  between  x. 
and  If  the  true  amplitude  probability  density  function,  f (x) ,  doesn't 

vary  "very  much"  in  this  interval,  then 


t 
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is  an  estimate  of: 


(7-2) 
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N(Xi+l  "  Xi} 


F(xm)  -  F(x.) 

(xi+l  -  V 
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If  the  successive  values  of  x(t)  are  independent  (that  is,  if  sufficient  time 
delay  has  occurred  between  samples),  the  probability  distribution  which 
describes  the  random  variable,  k.,  is  the  binomial  distribution.  If  p.  is 
the  true  probability  that  x^  <  x(t)  <x.+^,  then 

K  N-K. 

P(K  =  k  )  =  (  )  P.1  (1  -  P.)  1  (7-3) 

11  K. .  1  1 

1 

when  N  independent  samples  are  used  to  construct  the  histogram.  From  this  it 
may  be  shown  that 

E(K^)  =  p^N;  where  E  {•}  denotes  statistical  expectation 

so  that  E (P . )  =  p. 

l  l 


(7-4) 

(7-5) 


and  hence  the  estimate  of  the  probability  that  x(t)  falls  in  class  (x^, 
x^+j)  is  unbiased.  For  large  values  of  N  (in  the  present  case  N  will  always 
be  large  enough)  P(K^)  is  very  nearly  normal  and  may  be  used  to  put  a  confi¬ 
dence  interval  on  k^.  Recalling  from  the  theory  of  Binomial  Distributions 
that 

VAR(K  )  =  P.(1~P.)N  (7-6) 

it  can  be  said  with  8%  confidence  that  the  "true"  value  of  K.  is  contained  in 

l 

I - - 


If  p^  is  approximated  by  p^,  one  may  then  say  that  the  "true"  value  of  is 
contained  in 


h±Za\h  {1  '  T* 


where  Z  is  related  to  8  as  shown  in  Table  7-1. 
a 


Table  7-1.  Relationship  of  Z  to  6. 


.  percent 

Z  a/2 

90 

1.645 

95 

1.960 

98 

2.326 

99 

2.576 

The  equations  stated  in  this  section  may  be  used  to  estimate  the  number  of 
points  required  for  a  histogram  with  a  given  level  of  precision. 

If  the  fractional  error  for  the  i-th  class  is  defined  at  some  level  of 
confidence  8%  as 


e .  (8)  = 

i 


P  -  P 
i  i 


(7-9) 


then  using  the  result  of  Equations  (4)  and  (8)  it  may  be  shown  that 


ei(6) 


(7-10) 


where  p^  is  the  probability  that  a  given  x(t)  falls  in  the  interval  (x^, 
x^+f).  At  95%  confidence  one  can  then  generate  Table  7-2  which  shows  the 
total  number  of  data  points,  N,  required  to  obtain  a  given  error  as  a  function 
of  class  interval. 
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Table  7-2.  Data  Points  Required. 


Percent  Error 

20% 

10% 

5% 

2.5% 

.  2  a  \ 

width  /  p 

Class  Interval^  p  +  o 

Centered  1  p  +  2a 

at  '  p  +  3a 

1156 

1967 

9112 

111,222 

4624 

7868 

36443 

444,889 

18500 

31470 

145770 

1779556 

74000 

125883 

583084 

7182240 

Number  of  Points  N 

7.1.1  Guide  for  Obtaining  Histograms 

When  obtaining  histograms  one  should  observe  the  following  points: 

a.  For  the  histograms  to  be  useful  in  estimating  higher  order  moments 
(mean,  var,  etc.)  the  classes  must  be  sufficiently  fine.  Bendat(43) 
suggests  that  the  class  width  Ax  _<  .2o  everywhere  except  at  the 
tails  for  the  means  and  variances. 

b.  To  have  sufficient  data  for  testing,  mean  computation,  etc.  the 
experiment  should  be  planned  so  that  ki  >_  5  for  each  class.  To 
facilitate  this,  each  tail  may  be  made  a  class,  as  shown;  all 
other  intervals  have  widths  of  0.2o. 


c.  If  a  "histogram  analyzer"  is  used  with  a  countertype  laser  veloci- 
meter  processor  with  a  digital  output,  a  bias  error  may  occur  due 
to  the  mis-match  of  quantizing  steps  between  the  analyzer  and  the 
digital  device.  For  an  example  of  this  error,  see  Reference  44. 

d.  Use  enough  data  to  obtain  sufficient  accuracy  in  the  histogram, 
depending  on  the  purpose  for  which  this  data  will  be  used  where  2a/ 2 
is  related  to  R  as  shown  in  Table  7-1.  Tables  7-2,  7-3  and  7-4 
should  be  useful  in  determining  N  before  the  data  is  collected. 

GO 


7.2  STATISTICS  DERIVED  FROM  THE  HISTOGRAM 


Several  useful  statistics  may  be  derived  directly  from  the  histogram. 
Among  these  are  the  mean  (average  velocity)  and  the  moments  about  the  mean 
(central  moments).  The  set  of  central  moments  contains  the  variance  and 
standard  deviation  (turbulence  level),  skewness  and  "Kurtosis"  (peakedness). 

7.2.1  Mean 

Recall  that  the  mean  or  expected  value  of  a  discrete  random  variable  is 
defined  as: 

N 

Ux  =  E(X)  =  Z  x,  P(Xt  =  x4)  (7-14) 

i=l  1 


This  may  be  estimated  using  the  histogram.  Recall  that 


(7-15) 


is  an  estimate  of 


P(Xj.  _<  x(t)  <  *i+1),  the  probability  that  <x(t)  £  Xi+]/ 

Thus  if  the  class  intervals  are  small  enough,  one  may  estimate  the  mean  as 


Z 

all 

intervals 


x.,,  +  x.  k. 

(_i±l - -I)  _J_ 

v  2  N 


where  ( 


xi+l  +  Xi 


k. 

l 


and 


)  is  the  value  at  the  center  of  the  class  interval, 
is  the  number  of  counts  in  that  class, 
is  the  total  number  of  points  in  the  histogram. 


7.2.2  Moments  about  the  Mean 


The  £th  moment  about  the  mean  is  def ined  as 


Ef(x  -  uxr) 
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After  computing  x, 


one  may  estimate  the  itli  moment  about  the  mean,  m 


V 


as 


all  class 
intervals 


x.  .  .  +  x . 
.  l+l  1 


(7-16) 


These  estimators  are  biased,  but  for  large  n  (as  in  all  GE  LV  histograms) 

(45) 

this  bias  is  negligable 


7.2.3  Variance 

When  l  =  2,  Equation  7-16  gives  the  variance.  The  square  root  of  this 
estimates  the  standard  deviation  (or  turbulence). 


7.2.4  Skewness 

When  %  =  3,  one  obtains  the  third  moment  about  the  mean.  This  moment  is 
required  to  calculate  the  correlation  coefficient  between  jet  velocity  squared 
and  acoustic  pressure. 

If  m^  =  0,  there  is  no  skewness  (the  distribution  is  symmetric: 


If  m <  0,  the  distribution  skews  to  the  left: 


If  m-j  >  0  the  distribution  skews  to  the  right: 
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The  "Skewness"  (az  in  the  literature)  is  normalized  so  that  it  is  in¬ 
dependent  of  scale  and  is  calculated  as 


a3 


3/2 


(7-17) 


7.2.5  Kurtosis 

When  i  =  4,  the  moment  relates  to  the  peakedness  of  the  distribution. 

This  moment  is  required  to  calculate  the  correlation  coefficient  between  jet 

velocity  squared  and  acoustic  pressure.  The  normalized  version  of  this 

static  is  called  a,  and  is  calculated  as 
4 


_^4 
a4  2 
m2 


(7-18) 


If  the  distribution  is  Gaussian, 


a,  =  3.0 
A 


If  the  distribution  is  uniform 


1.8 


Sometimes  is  called  the  Kurtosis  and  other  times,  the  Kurtosis  is 
defined  so  that  the  Gaussian  distribution  has  zero  Kurtosis,  i.e.: 


=  <v3 


(7-19) 
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Both  appear  as  definitions  in  the  literature. 


7.3  ERROR  IN  MEAN  AND  TURBULENCE  MEASUREMENT 

The  statistical  error  in  the  mean  and  turbulence  measurements  obtained 
from  the  histogram  are  now  investigated.  If  the  histogram  is  approximately 
Gaussian  in  shape,  then  the  sample  mean  will  follow  Student's  t  distribution. 
If  one  lets  the  turbulence  level  (square  root  of  the  sample  variance)  be  s 
then  a  95%  confidence  statement  may  be  written  for  the  sample  mean  as: 

The  true  mean  is  contained  in 

x  ±  t  •  (S/flN) 

\  =  .025,  N 

(7-20) 

=  x  (1  ±  t  •  n/^b 

|  =  .025,  N 

where  n  =  s/x,  which  is  the  turbulence  ratio. 

If  the  fractional  error  is  defined  as: 

error  =  j x  -  u '  (7-21) 

u 

then,  with  95%  confidence,  it  may  be  said  that 

error  £  t  •  (i|/N)  (7-22) 

|  =  .025,  N 

For  convenience  the  error  is  tabulated  in  Table  7-3  as  a  function  of  N  and  n 
for  95%  confidence. 


(3d 


r 


Table  7-3.  Error  Percent  in  Mean  Measure¬ 
ment  with  95%  Confidence  as  a 
Function  of  N  and  n 


.2 

.1 

.05 

.025 

10 

14.1 

7 

3.5 

1.76 

20 

9.3 

4.7 

2.3 

1.2 

30 

7.4 

3.7 

1.9 

.93 

40 

6.3 

3.2 

1.6 

.8 

60 

5.0 

2.6 

1.3 

.65 

120 

3.6 

1.8 

.9 

.45 

In  a  similar  manner,  one  may  find  the  fractional  error  in  the  estimate 
of  the  turbulence,  S.  If  the  histogram  is  roughly  Gaussian, 


S2  (N  -  1)  2 

2  aN-1 


(7-23) 


2 

where  is  a  chi-square  random  variable  with  N-l  degrees  of  freedom.  Thus, 

a  confidence  interval  on  the  true  turbulence,  S,  is  then 


S 


.025,  N 


>  a  > 


.975,  N 


(7-24) 


with  95%  confidence.  The  worst  fractional  error  (since  the  tails  of  the  Chi- 
square  distribution  are  unsymmetric)  is  then 


error 


(N-l) 


(7-25) 


The  percent  error  at  95%  confidence  for  various  values  of  N  is  presented 


in  Table  7-4. 


Table  7-4.  Fractional  Error  in  Percent 
for  Turbulence  Estimate  as 
a  Function  of  N. 


N 

Error  Percent 

20 

31.5 

40 

21.8 

60 

17.8 

120 

12.6 

240 

9.12 

480 

6.45 

960 

4.56 

5000 

2.0 

25,000 

.89 

Examination  of  Tables  7-3  and  7-4  shows  that  considerably  more  data  must 
be  processed  for  a  turbulence  estimate  than  is  needed  for  a  mean  velocity 
estimate  of  equivalent  accuracy.  When  measuring  the  skewness  or  Kurtosis, 
similar  increases  in  data  for  equivalent  accuracy  will  be  required  over  the 
turbulence  estimate.  Also,  the  bias,  due  to  class  interval  width  will  in¬ 
crease  and  the  Bendat  rule  of  .2~>  will  probably  have  to  be  tightened.  It  is 
difficult  to  get  estimates  of  skewness  and  Kurtosis  repeatable  to  less  than 
several  percent. 

7 . 4  CORRELATION  FUNCTION  MEASUREMENT  WITH  THE  LASER  VELOCIMlTER 

The  measurement  of  correlation  functions  with  the  velocimeter  is  impor¬ 
tant  because  it  extends  the  usefulness  of  the  tool  allowing  it  to  perform 
measurements  of  turbulence  spectra  (obtained  by  Fourier  Transformation  of  the 
velocity  auto  correlation  function),  turbulence  eddy  convective  Mach  number, 
turbulence  length  scale,  the  spectra  of  turbulence  in  the  moving  frame  (this 
information  can  be  obtained  from  cross  correlation  between  the  velocity  at 
two  points  in  the  jet),  and  causality  (the  Lee-Ribner  approach  requires  the 
cross  correlation  between  the  jet  velocity  squared  and  the  far  field  pressure 
These  measurements  require  that  a  special  analysis  procedure  be  developed  for 


the  laser  velocimeter  that  accounts  for  the  random  sampling.  This  section 
presents  a  method  developed  at  General  Electric  that  has  been  successfully 
used  to  determine  the  important  aerodynamic  and  acoustic  parameters  enumerated 
above. 

The  procedure  is  developed  in  detail  for  the  autocorrelation  of  the  jet 
velocity,  x(t).  An  explanation  is  then  presented  to  show  how  the  estimator 
can  be  modified  for  determining  the  cross  correlation  of  the  velocities  at 
two  points  in  the  jet  and  for  determining  the  correlation  function  between 
the  jet  velocity  and  a  continuous  signal  such  as  that  obtained  at  the  output 
of  a  microphone. 

The  approach  considers  the  assumptions  of  wide  sense  stationarity  and 
ergodicity  made  in  all  methods  of  spectrum  analysis.  Wide  sense  stationarity 
prescribes  that  the  values  of  the  autocorrelation  function  (ACF)  are  functions 
of  the  time  between  the  signal  samples  and  not  the  times  the  samples  are 
taken.  Hence,  the  ACF  of  the  signal,  x(t),  is  written  as 


rxx(t)  =  E  {x(+)x(t+i)} 


(7-26) 


where  E { * }  denotes  expectation.  If  the  signal  x(t)  is  ergodic,  one  can 
interpret  the  expectation  as  an  average  over  time.  Therefore  an  acceptable 
estimate  of  R  (r)  for  some  particular  value  of  t  might  be 


Rxx(T)  X(tk+  T) 


(7-27) 


where  the  t^  are  M  random  times  at  which  samples  are  taken. 

Note  that  Equation  7-27  could  be  obtained  by  considering  the  relationship 
between  the  ACF's  of  the  signal  before  and  after  sampling^^’  .  If  the 
sampling  process  is  modeled  by  a  signal,  x(t),  consisting  of  impulses  occurring 
at  the  sampling  times,  then  the  sampled  version  of  x(t),  which  is  called  y(t) 
can  be  written 


where 


S(t)  =  ^  ^  u  (t-t  ) 
o  1 

i=-“> 


(7-29) 


and  t^  are  the  random  sampling  times.  If  it  is  assumed  that  s(t)  and  x(t) 
are  independent  random  processes. 


R 

yy 


(t)  +  R 


xx 


(t)  R„ 


(O 


(7-30) 


Then  the  ACF  of  x(t)  may  be  estimated  by  taking  the  ratio  of  estimates 
of  the  ACF's  of  the  sampled  version  of  the  signal,  y(t),  and  of  the  sampling 
signal,  s(t).  It  has  been  shown  (A3, 44)  that  the  ACF  of  s(t)  is  directly 
related  to  the  inter-arrival  probability  density  functions  of  the  sampling 
process.  These  are  determined  by  measuring  the  time  intervals  between  sample 
Substitution  of  the  traditional  estimators  into  the  ratio  gives  Equation  7-27 
The  derivation  of  Equation  7-27  sketched  above  requires  that 


R  (t)  i  0 
ss 


(7-31) 


for  all  t  of  interest.  This  requirement  is  equivalent  to  stating  that  there 
must  be  a  finite  probability  that  samples  are  taken  t  time  units  apart; 
otherwise  no  products  can  be  accumulated.  This  derivation  also  requires  that 
s(t)  and  y(t)  be  stationary  processes;  a  condition  which  is  not  really 
necessary . 

Finally,  apart  from  the  requirement  stated  in  Equation  7-31,  no  assump¬ 
tions  are  made  on  the  statistics  of  the  sampling  process.  The  ability  to 
construct  the  ACF  of  the  sampled  process  is  not  affected  by  the  mean  sampling 
rate.  There  is  no  "Nyquist  Criterion"  as  exists  in  traditional  sampled  data 
spectrum  analysis. 

ft  is  shown  below  that  the  ACF  reconstruction  technique  proposed  in 
Equation  7-27  yields  an  acceptable  estimate  of  RXX(T). 

In  particular,  it  will  be  shown  that  it  has  the  following  properties: 
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1 


It  Is  unbiased. 


2.  It  has  finite  variance. 

3.  The  variance  goes  to  zero  as  the  number  of  samples  in  the  estimate 
goes  to  infinity,  i.e.,  the  estimate  converges  in  the  mean. 

Thus  the  estimate  satisfies  the  conditions  necessary  to  be  called 
consistent. 


It  is  easy  to  show  that  R  (t)  will  be  an  unbiased  estimate  of  R  (t) 

xx  xx 


since 


M 

E(“xx('»  -5  Z  ■  “xx'1* 

k=l 


(7-32) 


One  may  now  derive  the  variance  of  the  ACF  estimate  given  in  Equation 
7-27  by  first  writing: 


M  M 


E{R  (t)} 

XX 


=  M2  2  2 

i=l  j-1 


E{X(t.)X(tj)X(t.+  1)x(tj+T)} 


■R2  Z  E(*2<ti)x2<tj>* 


i=l 


M  M 


1=1  .1=1 


E{X(t.)X(t.)X(t.+  i)X(t.+  T)} 
i  3  i  J 


(7-33) 


where  the  expectation 

t  ,+T  . 

J 


is  over  x(t);  and  t^  and  tj  given  samples  at  t^+i  and 


If  it  is  assumed  that  x(t)  is  jointly  normal  and  zero  mean,  the  result  for 
the  expected  product  of  four  normal  variables  can  be  used  to  simplify  the 
summations.  With 
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This  is  a  general  result  for  the  variance  of  the  estimator.  Written  in 
this  form,  the  first  term  represents  its  variance  if  all  of  the  product  pairs 
in  the  estimate  are  independent.  The  second  term  represents  the  increase  in 
the  variance  due  to  redundant  information  in  the  product  pairs.  Therefore, 
if  the  samples  are  taken  far  enough  apart  so  that 

R  (t.-t.)  =  0  all  t.  ,t.  (7-36) 

xx  j  y  if  J 

the  variance  of  the  estimate  is 
R2  (0)+R2  (') 

VAR{R  (T)}  =  — - r,-X-X -  (7-37) 

xx  M 

This  resuLt  may  still  he  a  good  approximation  to  the  true  variance  if 
the  mean  arrival  rate  of  product  pairs  (generally  much  less  than  the  mean 
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sampling  rate)  is  less  than  the  reciprocal  of  the  length  of  the  correlation 
function  of  x(t). 

While  Equation  7-35  gives  a  generic  result  for  the  variance  of  the  AFC 
estimate,  it  is  of  little  use  in  performing  "engineering"  calculations  of  the 
error  in  laser  velocimeter  correlation  functions.  However,  experiments  at 
General  Electric  have  shown  that  the  particle  arrivals  closely  follow  a 
Poisson  probability  distribution.  Assuming  a  Poisson  distribution  allows 
considerable  simplification  of  Equation  7-35  and  provides  some  simple  rules 
for  estimating  errors  in  correlation  measurements. 

Obtaining  the  result  for  the  Poisson  sampling  case  is  simplified  by  the 
fact  that  the  subset  of  samples  at  times  t^,  with  a  sample  in  [t^+x-A/2,t^ 
+T+A/2],  is  also  a  Poisson  process.  The  mean  arrival  rate  for  product  pairs 
can  be  shown  to  be 

X’  =  X2A  (7-38) 


where  A  is  the  "time  grid  width"  representing  the  amount  of  deviation  allowed 
in  selection  of  products  lag  r  apart,  and  X  is  the  mean  rate  of  the  original 
Poisson  process.  The  reason  for  the  time  grid  width.  A,  is  that  the  prob¬ 
ability  of  taking  two  samples  exactly  t  apart  is  zero.  Thus  in  any  practical 
realization  of  Equation  7-27  one  must  lump  together  those  samples  which  occur 
T-|<T<T+y  apart.  Averaging  product  pairs  at  close  to  the  same  lag  value 
results  in  the  addition  of  a  "time"  grid  jitter  error  to  the  estimate. 

Fortunately  it  may  be  shown  that  this  error  may  be  made  negligibly  small  if  A 
.  (44) 

is  small  enough 

To  simplify  the  derivation  of  the  bound,  one  can  assume  with  no  loss  of 
generality  that  the  t^'s  in  Equation  7-35  are  ordered  so  that 


t.>t.  when 
i  J 


One  may  then  write  Equation  7-35  in  terms  of  the  inter-arrival  density  functions 
of  the  Poisson  process.  These  are 
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(7-39) 


r  (AX2)ktk-1e-(AX  ^  , 
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Thus  Equation  7-35  becomes 
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/a,2,  k-1  ,2. 

(AX  )p  .-(AX  )p 
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(k-1)! 


Since  the  integrals  in  the  summation  are  all  positive,  one  may  write  an 
bound  on  this  result  as 

VAR{ R  ( r) }  <  Rxx(0)+Rxx(t) 

XX  - M - 
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'•o 
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[Rxx(y)+Rxx(,J+T)  Rxx(lJ-T)J 
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(k-1)! 


By  extending  the  summation  to  infinity,  an  identity  is  obtained  for  A^Ap 
Employing  Schwarz's  Inequality,  it  can  be  shown  that 


ao 
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R  (vj+t)R  (u-i)' 

XX  XX 
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Therefore  the  upper  bond  may  be  written: 


VAR{ R  (x)}< 
xx 
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This  is  the  upper  bound  for  the  variance  of  the  ACF  estimate  in  Equation  7-27 
when  the  sampling  is  the  Poisson  process. 


Consider  the  significance  of  the  two  terms  in  Equation  7-43.  Note  that 

2 

the  first  term  is  independent  of  A  and  that  the  second  term  depends  on  A  . 

Thus  for  the  low  mean  sampling  rates,  the  first  term  dominates  and  the  variance 
is  the  same  as  for  the  case  where  all  product  pairs  in  the  sum  are  independent. 
Now  if  A  becomes  very  large  so  that  the  second  term  dominates, 


VAR{R2  (t)} 


2A2A 

II 


R2  (p)dp 


(7-44) 


If  T  is  defined  as  the  length  of  time  over  which  the  observations  are  taken, 

L  2 
then  T^  =  M/A  A  and 


VAKIR^,))  If- 

L 


(t)  dx 


(7-45) 


If  the  mean  arrival  rate  is  very  large,  successive  samples  are  highly  correlated 
and  the  variance  is  reduced  as  the  observation  time  increases.  In  either 
case  the  variance  goes  to  zero  as  the  amount  of  data  processed  increases  so 
that  the  estimate  is  always  consistant. 


Usually,  the  laser  velocimeter ' s  mean  particle  arrival  rate  is  less  than 
10,000  particles/second.  For  most  jet  flows,  it  can  be  shown  that  the  first 
term  of  the  variance  equation  dominates  the  second.  Thus  for  engineering 
purposes,  the  variance  of  the  AFC  estimate  can  be  approximated  by 

.  R2  (0)+R2  (x) 

VAR{ R  (t)}  *  -55- - 55 -  (7-46) 

xx  M 


This  can  be  used  to  construct  confidence  limits  on  the  estimate  if  R  (t)  is 

xx 

assumed  to  be  roughly  Gaussian.  With  this  assumption,  a  98%  confidence 
statement  can  be  written  as 


73 


R  (t)  -2 

XX 


.33^ 


R  (t)  +2.331 
xx  1 


1  '2 

~2 

R 

0)+R 

(r) 

1  XX 

XX 

! 

M 

(T) 

< 

l'2 

'2 

k 

(0)+K 

(0 

XX 

x;: 
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The  estimation  procedure  Riven  in  Equation  7-46  can  be  extended  so  that 
the  cross  correlation  between  the  outputs  of  two  laser  velocimeters  can  be 
determined.  This  is  done  by  simply  accumulating  cross  product  pairs: 

H 

V'>  ’  ff  2  x(tk)y(V'>  <7-48> 

k=l 


The  procedure  can  also  be  used  to  cross  correlate  the  output  of  a  laser 
velocimeter  with  a  continuous  signal.  In  this  case,  let  y  (t^+T)  represent 
the  continuous  signal  that  is  physically  sampled  (using  a  sample  and  hold 
circuit  and  an  A/I)  converter)  at  time  t^+i .  For  a  particular  value  of  t,  a 
delay  circuit  can  be  used  to  cause  the  occurrence  of  this  sample  when  a 
particle  arrives  at  the  probe  volume  at  t^. 

The  procedures  described  above  have  been  used  successfully  at  General 
Electric  to  determine  the  auto  correlation  function  at  a  single  point  in  a 
jet,  the  cross  correlation  function  between  two  points  in  a  jet  and  the 
cross-correlation  between  the  jet  velocity  and  the  jet  velocity  squared  and 
far  field  microphone  signal. 


7. 5  MEASUREMENT  OF  EDDY  CONVECTION  VELOCITY 

Cross-correlation  data  obtained  between  the  outputs  of  two  laser  veloci 
meters  can  be  used  to  determine  the  velocity  of  turbulent  eddies  in  a  jet. 

To  do  this,  a  set  of  several  correlation  functions  are  obtained  witli  the 
laser  velocimeter  probe  volumes  at  different  spacings,  f ..  For  the  i*-'1 
correlation  function,  the  lag  value  corresponding  to  the  peak  in  the  correla 
t  ion  function,  t.  is  determined.  Thus 
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7T 


(7-49) 


t .  is  the  value  where 
1 


R  (t)  is  maximum 
xy 


Physically,  t ^  is  interpreted  as  the  time  required  for  the  eddy  to  travel 
from  the  first  laser  velocimeter  to  the  second.  Since  the  eddy  has  traveled 
the  probe  volume  separation  distance  in  this  time  interval,  an  estimate  of 
the  eddy  convection  velocity  is 


V. 

l 


(7-50) 


An  improved  estimate  of  eddy  velocity  is  obtained  by  considering  a  series  of 
such  measurements.  If  the  points  (i^,  f;^)  were  plotted  for  a  set  of  N  such 
measurements  on  a  piece  of  graph  paper,  the  eddy  convection  velocity,  V  , 
would  be  slope  of  the  best  line  that  could  be  drawn  through  the  points.  Thus 
when  determining  from  a  set  of  correlation  functions,  the  standard  linear 
regression  formulas  can  be  employed: 


M 


1=1 


and 


(7-52) 


The  convective  mach  number  can  be  estimated  by  dividing  by  the  appropriate 
jet  velocity. 
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7.6  DETERMINATION  OF  TURBULENCE  LENGTH  SCALE 


Turbulence  length  scale  can  also  be  determined  from  a  set  of  velocity 
cross  correlation  measurements  obtained  at  different  separation  distances. 
The  procedure  used  is  given  in  llinze^^.  Defining  the  ensemble  of  correla¬ 
tion  functions  as 


R 

xy 


(Y-i> 


(7-53) 


where  r  is  the  time  lag  and  t, .  is  the  separation  distance,  one  can  form  the 
longitudinal  correlation  coefficient 

R  (0,C.) 

KY)  =  ^ - ~L- -  (7-54) 

1  /R  (0)R  (0) 

xx  yv 


by  evaluating  each  correlation  function  obtained  at  T=0.  The  integral  length 
scale  is  then  given  from 


L 

x 


f  (Odr 


(7-55) 


This  integral  is  approximated  by  fitting  a  curve  of  the  form 


f  (O  =  a(exp\— — 


MsOy.) 


(7-56) 


(49) 


The 


to  tiie  data  points  f  (r  .  )  using  the  Conjugate  Gradient  procedure. 

]  1 

function,  f  (f),  is  then  numerically  integrated  using  the  Simpson's  rule 


procedure.  Gaussian  Quadrature  integration  could  also  be  employed  to  di  recti'. 


obtain  L  from  the  f(f .). 

x  i 


7  .  7  DETERMINATION  OF  LACRANG1AN  TIME  CORRELATION 

A  scries  of  two-point  velocity  correlations  at  different  separation 
distances,  can  also  he  used  to  determine  the  time  correlation  function  of 

i 

an  edflv  in  the  moving  reference  frame.  Since  the  separation  distance,  '. , 


and  the  lag  in  the  moving  frame  t^  are  related  by  the  convective  eddy  velocity, 
it  is  true  that  the  peak  values  of  the  cross-correlation  functions  represent 
respective  values  of  the  moving  frame  eddy  Correlation  function.  With  R^(t) 
defined  as  this  LaGrangian  correlation  function, 

RL(ti^  -  MAX{Rxy 

where  t^  =  the  value  of  t  where  the  maximum  occurs. 

The  spectrum  of  the  eddy  in  the  moving  frame  can  be  determined  by 
taking  the  Fourier  transform  of  the  above. 
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8.0  DESCRIPTION  OF  GENERAL  ELECTRIC'S  COMPUTATIONAL  ALGORITHMS 


Two  computer  programs  have  been  designed  to  perform  the  collection  and 
correlation  of  the  test  data.  The  program  INJET/INJF.T  handles  data  from  two 
laser  velocimeters  (LV1  and  LV2);  and  as  the  name  implies,  it  is  used  to  corre¬ 
late  velocity  data  from  two  points  in  the  jet  stream.  The  second  program, 
INJET/FARFLD,  operates  on  data  from  a  single  velocimeter  (LV1)  and  data  from 
a  far  field  pressure  transducer  (P2). 

Both  programs  employ  a  data  collection  and  summation  module  written  in 
assembly  language  for  efficiency  with  regard  to  execution  time.  The  input 
data  are  quantized  to  32  levels,  and  a  technique  called  quantized  product 
factoring  is  used  to  defer  the  multiplication  process  until  a  point  in  the 
data  processing  when  real  time  computations  are  not  required. 

Prior  to  the  start  of  the  data  summation  portion  of  both  programs,  an 

estimate  is  made  of  the  mean  and  variance  of  each  input.  This  is  done  by 

reading  128  samples,  summing  each  input  and  summing  the  square  of  each  input 

in  double  precision  integer  form.  The  mean  is  used  as  a  reference  or  DC 

level  for  all  subsequent  data  inputs  and  calculations.  The  variance  is  found 
2  -2 

from  EV  /128  -  V  and  is  used  in  determining  the  cell  size  to  be  used  for 
quantizing.  A  cell  size  equal  to  3/16 /variance  assumes  that  32  levels  will 
encompass  99%  of  the  total  input  range  for  a  normal  distribution. 

Although  the  two  programs  are  similar  in  structure,  the  details  of  the 
data  collection  and  summation  differ.  For  clarity,  the  two  programs  will 
be  discussed  separately. 

8.1  QUANTIZED  PRODUCT  FACTORING 

Quantized  product  factoring  might  best  be  described  by  the  following: 

The  normal  expression  for  the  sum  of  products  is: 
m 

Z  xtyi 

i=l 
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If  the  variables  x  and  y  are  quantized  to  N  levels,  there  are  N  finite  values 
that  the  quantized  variables,  X  and  Y,  can  take  on.  This  means  that  the 
product  can  now  be  factored,  and  further  that  the  quantized  values  can  be 
used  as  an  index  or  pointer  into  a  vector.  The  quantized  factored  form  of 
the  sum  of  products  is 


Quantized  factorization  means  that  summations  may  be  used  instead  of  multi¬ 
plications  during  that  portion  of  the  data  processing  which  is  being  done  in 
real  time. 

8.2  IN JET/ IN JET  DATA  COLLECTION 

All  velocity  data  is  collected  using  an  interrupt  technique.  When  a 
velocity  data  point  is  obtained,  that  is,  the  LV  circuitry  has  Indicated  that 
the  A/D  conversion  is  complete,  an  interrupt  is  generated  in  the  computer. 

Four  words  are  read  from  the  digital  input  device  in  response  to  each  inter¬ 
rupt,  Tl,  VI,  T2  and  V2  in  that  order.  T1  and  T2  are  the  external  clock  times 

at  the  start  of  the  A/D  conversion,  VI  and  V2  are  the  A/D  output  values. 

Since  either  LV  can  generate  the  interrupt,  the  sign  of  the  velocity  word 
causing  the  interrupt  is  forced  to  a  1  to  indicate  valid  data. 

Input  data  for  each  LV  is  stored  on  stacks  of  1024  words  each.  This 
allows  room  for  512  time  and  velocity  pairs.  A  counter  keeps  track  of  the 
remaining  room  on  each  stack,  and  a  pointer  locates  the  next  open  location  on 
each  stack.  When  these  pointers  reach  the  bottom  of  the  stack,  they  are  reset 

to  the  top  of  the  stack  as  long  as  the  counter  indicates  space  available. 


The  time  value  is  placed  on  the  data  stack  directly.  The  velocity 
value  put  on  the  stack  is  a  quantized  velocity  relative  to  the  mean  or  DC 
level,  i.e. ,  Q1  =  (VI  -  V1)/CS1  where  CS1  is  the  cell  size  for  LV1. 

This  data  input  process  goes  on  independent  of  the  correlation  process 
which  will  free  locations  on  the  data  stacks  as  the  data  is  used,  with  the 
oldest  data  being  processed  first. 

8.3  VELOCITY  CORRELATION  AND  HISTOGRAM  ALGORITHMS 

The  correlation  process  is  a  pairing  of  each  data  point  on  one  stack  with 
all  data  points  on  the  other  stack  whose  time  differnce  (AT)  satisfies  some 
predefined  criteria.  A  total  of  512  lag  values  (AT5^-^)  was  chosen,  result¬ 
ing  in  a  range  for  AT  of  -256  <_  AT  <_  +255. 

A  rectangular  array,  32  columns  by  512  rows,  is  used  to  sum  the  results 
of  the  pairing  process.  The  AT  value  is  used  as  a  row  pointer,  since  each 
row  corresponds  to  a  specific  lag  value  from  -256  to  +255.  Each  column  cor¬ 
responds  to  a  quantizing  level  from  -16  to  +15.  When  a  pair  of  data  points 
is  found  which  meets  the  AT  criteria,  one  quantized  velocity  (say  Ql)  is 
used  as  a  column  pointer,  AT  as  the  row  pointer,  and  Q2  is  added  into  the 
specified  array  element.  This  is  using  the  quantized  factorization  technique 
described  previously.  Histograms  for  each  velocity  are  formed  by  using  the 
paired  Q  values  as  pointers  into  vectors  32  elements  long  and  incrementing 
the  corresponding  elements.  The  number  of  occurrences  of  each  lag  value  is 
generated  in  a  similar  manner  by  using  AT  as  a  pointer  into  a  512  element 
vector  and  incrementing  the  indicated  element. 

When  data  points  on  one  stack  do  not  pair  with  any  points  on  the  other 
stack,  the  locations  they  occupy  become  available  for  input  again.  This  pair¬ 
ing  and  summing  process  is  continued  until  some  predetermined  number  of  data 
points  has  been  processed,  usually  20,000  to  30,000  points;  then,  all  data 
input  and  processing  is  halted.  This  is  done  to  prevent  an  overflow  from 
occurring  in  the  summing  array. 


Another  interrupt  is  used  to  detect  the  overflow  of  the  external  clock. 
When  the  clock  overflows,  data  input  is  temporarily  halted  until  all  the  data 


currently  on  the  data  stacks  have  been  correlated,  then  the  data  input  and 
processing  is  resumed  again.  This  is  done  because  the  time  information  be¬ 
comes  ambiguous  after  a  clock  overflow  unless  a  double  precision  clock 
counter  is  used. 

8.4  DATA  COMPRESSION  AND  INTERMEDIATE  DISPLAYS 

When  the  required  number  of  data  points  has  been  examined  and  process¬ 
ing  halted,  the  summing  array  is  compressed  to  a  double  precision  vector. 

This  is  done  by  multiplying  each  column  of  the  array  by  its  corresponding 
quantized  velocity  value  and  summing  across  the  columns.  In  this  manner,  the 
more  time  consuming  multiplication  process  is  done  at  a  point  where  execution 
time  is  not  critical.  The  double  precision  vector  now  contains  the  sum  of 
the  quantized  product  pairs  for  each  lag  vlaue  (AT). 

New  estimates  of  the  mean,  variance  and  cell  size  are  calculated  for 
each  velocity  from  the  histogram  dat.'.  The  histogram  data  are  then  normalized 
and  displayed  on  a  scope  for  operator  monitoring. 

The  current  cross-correlation  contribution  is  calculated  from  the  double 
precision  sum  of  the  product  pairs  and  the  previous  mean  and  cell  sizes. 

The  total  run  correlation  coefficient  is  then  updated  along  with  the  total 
mean  and  variance  values.  A  normalized  correlation  coefficient  is  then  cal¬ 
culated  and  displayed  on  the  scope.  The  number  of  data  points  in  each  histo¬ 
gram  along  with  the  total  number  of  data  points  processed  in  the  run  is  also 
displayed  on  the  scope. 

A  new  DC  level  is  calculated  using  the  latest  value  of  the  mean,  all 
buffers  are  cleared  and  the  data  collection  and  summing  process  is  resumed 
again. 

This  processing  of  data,  in  groups  of  20,000  to  30,000  points,  is  re¬ 
peated  until  a  total  of  approximately  300,000  points  have  been  sampled,  and 
then  the  run  is  terminated.  Although  the  segmented  processing  is  required  to 
prevent  overflows  in  the  summing  buffer,  it  also  allows  the  quality  of  the 
data  collection  to  be  monitored  at  intervals  during  the  run.  This  means  that 
bad  data  can  be  detected  earlier  and  the  run  aborted,  rather  than  at  the  end 
of  a  complete  run. 
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8.5  INJET/FAR FLD  DATA  COLLECTION 


Pressure  data  are  obtained  on  an  interrupt  basis  from  a  separate  A/D 
converter  that  is  part  of  the  11/45  computer  system.  The  conversions  are 
initiated  by  signals  from  the  external  clock  that  is  part  of  the  laser  veloc- 
imeter  circuitry.  When  the  conversions  are  complete,  an  A/D  interrupt  is 
generated.  The  external  clock  rate  is  divided  by  4  so  that  the  pressure 
interrupts  will  occur  at  a  20  KHz  rate.  The  11/45  clock  is  set  up  to  gen¬ 
erate  a  clock  interrupt  when  the  external  clock  overflows.  A  software  clock 
is  used  internally  in  the  program  and  is  incremented  every  time  a  pressure 
interrupt  is  serviced. 

Velocity  data  are  obtained  through  the  same  digital  input  device  that  is 
used  for  the  INJET/INJET  program,  except  that  only  one  velocimeter  is  used. 

A  device  flag  is  set  rather  than  generating  an  interrupt  when  a  new  sample  of 
velocity  and  time  is  ready. 

The  down  stream  separation  between  the  velocimeter  and  the  microphone  is 
converted  to  a  delay  time  in  clock  count  units  using  the  speed  of  sound  in 
the  test  cell  and  the  clock  frequency.  This  value  is  reduced  by  256  counts 
so  that  a  time  window  for  correlation  is  centered  at  the  mircophone. 

When  each  velocity  and  time  pair  is  read,  the  velocity  is  placed  on  a 
delay  stack,  and  the  time  value,  increased  by  the  delay  counts,  is  placed  on 
a  start  time  stack.  A  delay  stack  counter  is  decremented  to  show  one  less 
space  is  available  on  the  delay  stack.  When  the  software  clock  matches  the 
start  time,  the  corresponding  velocity  value  is  removed  from  the  delay  stack, 
512  counts  are  added  to  the  start  time  to  form  a  stop  time,  and  this  value  is 
placed  on  a  stop  time  stack.  The  velocity  is  then  paired  with  the  next  512 
pressure  values  as  they  are  read  in  the  pressure  interrupt  routing. 

8.6  VELOCITY-PESSURE  CORRELATION  AND  HISTORGRAM  GENERATION 

The  correlation  of  each  velocity-pressure  pair  is  done  entirely  in  the 
pressure  interrupt  routine.  The  quantized  factoring  technique  is  utilzed  to 
defer  the  actual  multiplication  operation.  A  quantized  velocity,  Ql,  is 
formed  as  indicated  in  Section  8.2  and  is  used  as  a  column  pointer  into  the 
summing  array. 
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The  summation  is  done  entirely  within  the  pressure  interrupt  routine. 

In  response  to  each  pressure  interrupt,  the  A/D  value  is  read  and  a  quantized 
pressure,  Q2,  is  formed  in  a  manner  similar  to  the  velocity  quantizing.  In 
this  case,  no  subtraction  of  a  DC  level  is  made;  since,  by  definition,  the  DC 
pressure  level  is  zero. 

Implied  in  the  definition  of  the  time  window  and,  therefore,  the  start  and 
stop  times  is  that,  at  the  start  time  At  =  -256  clock  counts  and  at  the  stop 
time  At  =  +255  clock  counts.  Therefore,  when  a  velocity  value  is  "started". 

At  can  be  thought  of  as  pointing  to  the  first  row  in  the  summing  array.  The 
quantized  pressure  Q2  can,  therefore,  be  summed  into  the  first  element  of  the 
column  pointed  to  by  the  quantized  velocity  Ql. 

Since  each  pressure  interrupt  will  occur  one  clock  count  later,  the  next 
Q2  value  will  be  summed  into  the  next  element  of  the  same  column  which  corre¬ 
sponds  to  At  =  -255.  This  summing  process  is  continued  until  the  clock  time 
matches  the  stop  time  for  this  velocity  value  or  512  clock  counts  later. 

The  pressure  interrupt  routine  is  coded  to  handle  four  of  these  summing  and 
auto-indexing  sequences,  which  are  referred  to  as  accumulators.  Therefore, 
four  different  velocity  values  can  be  paired  with  the  same  pressure  value. 
Twelve  velocities  may  be  waiting  on  the  delay  stack. 

When  a  velocity  is  "started",  it  is  removed  from  the  delay  stack  and  the 
delay  stack  counter  is  incremented  to  indicate  an  available  location  for  a 
new  velocity.  A  check  is  made  to  see  if  an  accumulator  is  available.  If  all 
accumulators  are  in  use,  the  velocity  value  is  discarded.  However,  if  an 
accumulator  is  available,  the  quantized  velocity,  Ql,  is  used  to  increment 
the  velocity  histogram  as  described  in  Section  8.3,  the  stop  time  is  calcu¬ 
lated  and  stored  on  the  stop  stack,  and  the  processing  counter  is  incremented. 
This  is  the  counter  that  is  used  to  halt  processing  when  20,000  or  30,000 
points  have  been  processed.  The  accumulator  chosen  is  loaded  with  the  Ql 
value  and  enabled  so  that  during  the  next  512  pressure  interrupts  it  will  be 
active. 

When  the  clock  value  matches  a  stop  time,  the  corresponding  accumulator 
is  turned  off  and  the  last  quantized  pressure  value  is  used  to  increment  the 
pressure  histogram. 
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External  clock  overflows  generate  an  interrupt  which  causes  the  software 
clock  to  be  reset  to  zero.  To  prevent  this  from  happening  while  any  accum¬ 
ulator  is  running,  a  velocity  value  will  never  be  put  on  a  delay  stack  if  the 
delay  time  plus  the  accumulated  time  will  exceed  the  remaining  time  in  the 
current  clock  cycle. 

8.7  DATA  COMPRESSION  AND  DISPLAY 

The  compression  of  the  summing  array  to  a  double  precision  vector  con¬ 
taining  £Q1*ZQ2  is  accomplished  as  described  in  Section  8.4  for  the  INJET/ 

2 

INJET  program.  An  additional  step  is  introduced  so  that  the  sum  of  Q1  *Q2  is 
formed  for  each  lag  value. 

The  remainder  of  the  calculations  and  display  generations  is  essen¬ 
tially  identical  to  the  procedure  described  for  the  INJET/INJET  program. 


9.0  EXPERIMENTAL  ARRANGEMENTS  USED  FOR  LASER  VELOCIMETER  MEASUREMENTS  - 
GENERAL  ELECTRIC  ANECHOIC  JET  NOISE  FACILITY 


This  section  contains  a  brief  sketch  of  the  physical  experimental 
arrangement  of  General  Electric  Laser  Velocimeter  system  as  configured  for 
the  Anechoic  Jet  Noise  Test  Facility. 

The  nozzle  under  test  and  two  LV  systems  are  similar  to  what  is  shown  in 
Figures  9-1  and  9-2.  The  laser  beams  are  highlighted  above  the  nozzle  by 
nature  of  the  seeding.  In  the  right  of  Figure  9-1  the  two-velocimeter  pack¬ 
age  may  be  seen  mounted  on  its  main-three- axis  (X-Y-Z)  actuator.  The  LV  to 
the  right  is  installed  on  an  independently  operated  sub-three-axis  (X'-Y'-Z') 
actuator  that  controls  the  spacing  of  the  twin  LV  measurement  control  col¬ 
umns.  Control  and  data  acquisition  consoles  are  shown  in  the  photographs  of 
Figures  9-3  and  9-4.  The  minicomputer  is  in  the  three  low  cabinets  to  the 
right  in  Figure  9-3.  The  keyboard  input  and  printer  is  at  the  far  right. 

The  Mark  II  signal  processor  is  in  the  left-hand  cabinet,  and  the  Mark  III 
processor  is  on  the  lower  shelf  of  the  cart.  This  photograph  was  taken 
before  equipment  locations  were  finalized.  The  console  of  Figure  9-4  is 
mainly  for  nozzle  operation,  with  the  exception  of  the  right-hand  cabinet, 
that  houses  the  LV  actuator  controls. 
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10.0  RESULTS  OF  MEASUREMENTS 


The  General  Electric  Laser  Velocimeter  system  has  undergone  a  series  of 
electronic  and  automated  data  processing  developments  over  the  last  several 
years.  For  example,  within  the  last  two  years,  new  improvements  in  signal 
processing  and  data  processing  techniques  have  increased  data  taking  capability 
by  a  factor  of  90  for  two-point,  space-time  type  of  measurements.  Sections 
3.0  through  8.0  describe  many  of  the  existing  improvements  and  techniques 
necessary  for  performing  advanced  real-time  measurements.  References  30,  31, 

50  and  51  provide  some  of  the  reference  material  in  which  has  been  demon¬ 
strated  unique  LV  diagnostic  capability  for  heated  high  velocity  exhaust 
jets. 

The  first  item  discussed  in  this  section  is  a  summary  of  some  highlights 
of  previously  acquired  mean  velocity  and  turbulent  velocity  measurements 
performed  on  heated  supersonic  jets  and  of  turbulence  spectra  measurements. 
Secondly,  results  are  presented  of  more  recently  obtained  two-point, 
space-time,  in-jet  measurements;  and  in-jet  to  far-field  cross-correlation 
measurements  performed  on  a  conic  nozzle  and  a  coannular  plug  nozzle. 


10.1  MODEL  SCALE  MEAN  VELOCITY  AND  TURBULENT  VELOCITY  MEASUREMENTS 

The  laser  velocimeter  measurements  reported  herein  were  made  in  General 
Electric's  Jet  Noise  Acoustic  Test  Facility.  Two  nozzles  were  used  in  the 
tests:  a  converging-diverging  nozzle  producing  shock-free  supersonic  flow, 
and  a  conical  nozzle  used  for  subsonic  tests  and  shock  structure  measure¬ 
ments.  The  exit  diameter  of  both  nozzles  was  6  inches. 


The  velocity  measurements  were  obtained  as  probability  density  distri¬ 
butions  and  the  data  were  reduced  by  using  statistical  techniques  discussed 
in  Sections  7.0  and  8.0  to  compute  the  mean  values  and  standard  deviations  (mean 
velocity  and  turbulent  velocity),  respectively.  Typical  mean  and  turbulent 
laser  velocimeter  (LV)  measurements  for  an  ambient  subsonic  jet  are  compared 
with  hot  film/hot  wire  measurements  in  Figures  10-1,  a  and  b.  LV  measurements 
of  mean  velocity  and  turbulent  velocity  were  then  obtained  for  the  high- 
temperature  (1500°  R) ,  fully-expanded  supersonic  (Mj  =  1.5)  convergent-divergent 


90 


Instrument 


•  O  Laser  Velocimeter 
O  Hot  Film 
Hot  Wire 

♦O  Laser  Velocimeter 


Ambient 
Ambient 
Ambient 
1500°  R 


Solid  Symbols:  r/D 
Open  Symbols:  r/D 


p  0.6 


0.2  l-a.  Mean  Velocity  -  Ambient 
I  Subsonic  Jet.  I  I  I 


c.  Mean  Velocity  -  Hot,  Super 
.  sonic.  Shock-Free  Jet.  , 


0.20  |-b.  Turbulent  Velocity  -  Ambient 
I  Subsonic  Jet.  .  ■  i  I  I 


nozzle.  The  results  of  these  measurements  are  presented  in  Figures  10-lc  and 
d.  In  Figures  10-la  and  c,  the  existence  of  the  entended  potential  core  for 
the  supersonic  jet  is  evident.  The  differences  in  turbulence  distributions 
between  the  subsonic  and  supersonic  jet  can  be  observed  by  comparing  Figures 
10-lb  and  d. 

Further  experimental  measurements  were  obtained  using  the  converging- 
diverging  shock-free  nozzle  flow  and  the  shocked  conic  nozzle  flow.  For  both 
cases,  the  exit  Mach  number  of  the  jet  was  1.55  and  the  total  temperature  was 
1500°  R.  Radial  profiles  of  mean  axial  velocities  were  obtained  for  various 
axial  locations,  X,  in  the  jets  of  both  nozzles.  The  results,  presented  in 
Figures  10-2a  and  b,  are  plotted  using  the  parameter,  n,  which  collapses  all 
of  the  data  onto  one  curve  and  clearly  defines  the  jet  boundary  and  the  shear 
layer  at  r  =  r0.  For  the  conic  nozzle  case.  Figure  10-2b,  the  under-expansion 
in  the  jet  boundary  is  apparent  as  the  value  of  the  normalized  mean  velocity 
becomes  greater  than  1.0  in  the  expansion  region  near  the  jet  exit  plane  (X/D 
=  2). 

Radial  profiles  of  normalized  turbulence  velocity  for  the  converging- 
diverging  nozzle  are  shown  in  Figure  10~2c  as  a  function  of  the  same  variable, 
n,  (used  in  the  previous  figures)  and  was  again  used  here  to  collapse  all  the 
profiles  onto  one  curve.  The  peak  turbulence  (maximum  shear)  is  clearly 
observed  to  occur  at  the  jet  boundary,  r  =  rG. 

To  complement  the  above  conic  nozzle  measurements,  a  series  of  LV  measure¬ 
ments  was  also  made  to  define  the  conic  nozzle  shock  structure.  Figure  10-3 
shows  the  measured  axial  mean  velocity  and  turbulent  velocity  distributions 
along  the  centerline  of  the  jet.  The  expansion  and  contraction  of  the  axial 
velocity  in  the  shock-flow  region  are  clearly  shown.  For  reference  purposes, 
predicted  velocity  is  also  shown.  The  amplification  of  turbulent  velocity 
due  to  the  shock  is  also  evident. 

10.2  TURBULENCE  SPECTRA 

To  construct  turbulence  spectra  from  the  discontinuous  output  of  the  LV 
presents  a  difficulty  not  encountered  with  continous  type  instrumentation. 
Conventional  point  spectra  estimate  techniques  assume  that  all  values  of  the 
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Mean  Velocity;  C/D  Nozzle,  b.  Mean  Velocity;  Conical  Nozzle, 

Shock-Free  Data,  Mj  =  1.55,  Shocked  Data,  Mj  =  1.55, 

Tt  =  1500°  R.  TT  =  1500°  R. 


T!  =  (r-r0)/(X-D) 


=  (r-r0)/(X-D) 


c.  Turbulence  Velocity;  C/D  Nozzle, 
Shock-Free  Data,  Mj  =  1.55, 

Tt  =  1500°  R.  J 


n  =  (r-r0)/(X-D) 


Figure  10-2.  LV  Measured  Radial  Profiles  of  Normalized  Mean  Velocity  and 
Normalized  Turbulent  Velocity. 


Velocity,  fps 


1 


input  signal  are  known  in  the  analysis  interval.  Such  knowledge  is  not 
available  at  the  LV  output.  To  overcome  this  problem,  it  was  necessary  to 
construct  the  autocorrelation  function  of  the  velocity  signal,  then  obtain 
the  spectrum  as  its  Fourier  transform.  Sections  7.0  and  8.0  describe  the 
theoretical  and  computational  considerations  for  these  types  of  measurements. 
Figure  10-4a  shows  a  comparison  of  the  LV  reconstructed  spectra  with  that 
using  a  hot  film  for  a  subsonic  jet.  Figure  10-4b  shows  the  LV  measured 
turbulent  spectra  for  a  sonic  (Mj  •  1.0)  heated  (total  temperature  of  1500°  R) 
jet.  These  results  were  quite  encouraging  in  regard  to  the  real-time 
measurement  studies.  These  results  gave  more  confidence  to  pursue  cross¬ 
correlation  type  real-time  studies. 

10.3  TWO-POINT,  SPACE-TIME,  IN-JET,  CROSS-CORRELATION  MEASUREMENTS 

In  order  to  develop  the  technology  of  using  the  laser  velocimeter  for 
performing  advanced  type  turbulence  measurements,  a  combination  of  signal 
processing  development,  theoretical  analysis,  computer  software  data  process¬ 
ing  techniques  and  laser  velocimeter  actuation  equipment  development  was 
necessary.  Sections  3.0  through  9.0  describe  the  extent  of  development 
necessary.  The  measurement  of  two-point  cross  correlations  with  a  non¬ 
contact-type  measuring  device  such  as  a  laser  velocimeter  has  never  been 
performed  before.  The  full  accomplishment  of  such  techniques  will  enable  the 
measurement  and  location  for  an  exhaust  jet's  aerodynamic  noise  source  and 
will  eventually  lead  to  a  new  diagnostic  device  for  evaluating  and  better 
designing  exhaust  nozzle  suppression  devices. 

Discussed  below  are  results  of  a  series  of  two-point,  space-time,  in-jet 
type  measurements  in  a  conic  nozzle  and  a  coannular  plug  nozzle.  The 
resultant  measurements  could  be  used  in  the  aero-acoustic  indirect  models 
described  in  Section  6.0. 

10.3.1  Conic  Nozz  l_e  Me  as  u  rements 

Although  all  the  computational  schemes  and  algorithms  for  performing  the 
two-point,  space-time,  cross-correlation  measurements  are  covered  in  detail  in 
Sections  7.0  and  8.0,  it  is  instructive  to  briefly  review  the  concepts  here. 
Sketches  shown  in  Figures  10-5a  and  b  serve  the  purpose  for  the  following 


Comparison  of  Laser  Velocimeter  and 
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Figure  10-4.  LV  Measured  Axial  Turbulent  Velocity  Spectra  for  Ambient  and  Heated  Jets 
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b.  Turbulent  Properties  Defined  by  In-Jet  Two-Part 
Space-Correlat ion 


Figure  10-5.  Sketches  Illustrating  the  Coordinate  System  and  Turbulent 
Properties  Defined  by  Two-Point  Space-Time  Cross 
Correlation  Measurements. 


discussion  of  results.  Shown  in  figure  10- 5b  .ire  the  concepts  necessary  to 
compute  the  turbulent  length  scale,  Lx,  the  convection  velocity,  Vc,  and  the 
spectra  of  the  moving  eddy  turbulence.  Measurements  of  the  turbulent  structure 
properties  were  performed  in  General  Electric  Anechoic  Jet  Noise  Test  Facility 
described  in  Section  9.0.  Figures  9-1  and  9-2  are  photographs  of  the  two- 
laser  velocimeter  system  built  for  the  following  measurements. 

For  the  conic  nozzle  measurements,  two  nozzle  conditions  were  tested: 
Subsonic  (Mj  =  0.5)  ambient,  and  sonic  (Mj  =  1.0)  heated  (T’p  =  1600°  R) . 

Figure  10-6  shows  a  sample  of  six  two-point,  space-time,  cross  correlations  for 
the  subsonic  ambient  conic  nozzle.  The  axial  spacings  for  this  test  were 
f,/D  =  0,  0.25,  0.45,  0.6,  0.85  and  1.0;  where  f,  is  the  axial  separation 
distance  between  the  two  LV  control  volumes,  and  D  is  the  diameter  of  the 
conic  nozzle  (D  =  4.64  inches).  The  ordinate  on  each  scale  is  the  normalized 
cross-correlation  function,  Rxt>  while  the  abscissa  is  the  normalized  time 
scale,  t  Vj/D  (where  x  is  the  time,  Vj  the  ideal  exhaust  velocity).  From 
such  a  series  of  cross-correlation  measurements,  the  convective  speed  and 
length  scale  can  be  determined  as  illustrated  in  Figure  10-5b.  Table  10-1 
summarizes  the  results  of  the  conic  nozzle  measurements. 

To  illustrate  the  computation  of  the  convection  velocity,  Vc,  Figure  10- 
7  was  constructed  from  the  test  results  shown  on  Figure  10-6.  These  results 
show  that,  at  an  axial  location  of  4  diameters  at  x/D  =  0.25,  the  convection 
velocity,  Vc,  is  0.723  Vj  for  both  the  ambient  subsonic  jet  and  the  sonic 
treated  jet.  Summarized  in  Figures  10-8a  and  b  is  the  LV  measured  axial 
variation  of  the  turbulent  length  scale  and  the  radial  variation  of  convec¬ 
tion  velocity  for  the  measurements  taken.  The  test  results  shown  in  Figure 
10-8a  show  that  the  turbulent  length  scales  for  the  subsonic  unheated  and  the 
sonic  heated  jet  are  much  the  same.  There  is,  however,  a  noticeable  radial 
variation  of  length  scale,  but  it  is  not  large.  Also,  the  classically  quoted 
value,  Lx/D  =  0.12  X/D,  is  at  variance  with  the  above-measured  length 
scales,  as  well  as  some  University  of  Toronto  hot  wire  measurements.  Instead 
of  Lx/D  =  0. 12  X/D,  perhaps  a  more  representative  equation  would  be  Lx/D 
-  0.03  X/D. 

The  radial  variation  of  convection  velocity  for  the  subsonic  unheated 
and  the  sonic  heated  jet  measurements  are  summarized  in  Figure  10-8b.  These 
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Table  10-1.  Summary  of  Turbulent  Properties  Measured  on  a 
Conical  Nozzle. 


Subsonic  Ambient  Conical  Nozzle  Test  Results  (Mj  =  0.5;  Total 
Temperature  =  519°  R) 


X/D 

r/D 

V,  fps 

u'/v 

Vc/Vj 

Lx/D 

Lx ,  in .  1 

4 

0 

531.4 

0.043 

0.75 

0.057 

0.26 

0.25 

516.3 

0.11 

0.  723 

0.091 

0.42 

0.5 

152 

0.471 

0.414 

0.172 

0.8 

8 

0 

485.3 

0.132 

0.809 

0.211 

0.98 

0.25 

390.0 

0.215 

0.580 

0.223 

1.03 

14 

0 

304.9 

0.195 

0.487 

0.396 

1.84 

b.  Sonic  Heated  Conical  Nozzle  Test  Results  (Mj  =  1.0;  Total 
Temperature  =  1600°  R) 


X/D 

r/D 

V 

U’  /V 

Vc/Vj 

Lx/D 

Lx,  in. 

■ 

0 

0.25 

1739.7 

1511.2 

0.048 

0.10 

0.74 

0.797 

0.0428 

0.0934 

0.2 

0.43 

where:  X  =  Axial  location  measured  from  the  jet  centerline 

r  =  Radial  location  measured  from  the  jet  centerline 
Vc  =  Measured  turbulent  convective  speed 
V  =  Measured  local  mean  velocity 
Vj  =  Ideally  expanded  jet  velocity 
U'  =  Measured  axial  turbulent  velocity 
U'/V  =  Local  turbulence  intensity 
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Figure  10-6.  LV  Measured  Two-Point  Cross  Correlation  for  a  Range  of  Axial  Separation 
Distances  for  a  Conical  Nozzle. 


Vj/o 


Figure  10-7.  Illustration  of  Method  for  Determining 
Convective  Results  of  the  Two-Point 
Space-Time  Cross  Correlations. 
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b.  Radial  Variation  of  Convective  Speed 
for  a  Conical  Nozzle  as  Measured  by 
a  Laser  Velocimeter. 


Figure  10-8.  LV  Turbulent  Structure  Measurements  for  a  Conical  Nozzle 
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results  show  that  unheated  and  heated  convective  speeds  are  similar  in  value, 
but  there  does  exist  an  axial  and  radial  variation  in  value  throughout 
the  exhaust  jet. 

10.3.2  Coanmil_ar  PI ug_  Nozzle  Test  Rosul  ts* ** 

A  second  series  of  two-point,  space-time  measurements  were  performed  on  a 
second  nozzle  to  explore  the  type  of  phenomenon  and  problems  that  are  encoun¬ 
tered  with  nozzle  systems  that  are  different  than  the  simple  conic  nozzle. 

The  second  nozzle  selected  was  a  coannular  plug  nozzle  configuration  illus¬ 
trated  in  Figure  10-9.  Also  shown  on  Figure  10-9  are  the  two  nozzle  test 
conditions  for  which  the  cross-correlation  measurements  were  performed.  The 
coannular  plug  nozzle  tested  was  a  high  radius  ratio  (R°  =  0.85),  low  area 
ratio  (A^  =  0.33)  nozzle  tested  acoustically  under  Contract  NAS3-19777.  This 
particular  nozzle  has  been  found  to  exhibit  favorable  jet  acoustic  reduction 
characteristics  relative  to  a  conic  nozzle  at  an  equivalent  specific  thrust. 
The  test  conditions  chosen  were  such  that,  for  the  first  test  point,  the  outer 
and  inner  streams  were  slightly  less  than  critical;  for  the  second  test 
point,  the  inner  stream  was  maintained  at  the  same  conditions  as  for  test 
point  1,  but  the  outer  stream  was  raised  to  a  super  critical  pressure  ratio. 
Table  10-2  summarizes  the  results  of  all  the  two-point,  space-time  measure¬ 
ments  taken  for  this  test  series. 

For  this  coannular  plug  nozzle,  in-jet,  cross-correlations  were  generally 
taken  at  three  radial  locations  for  each  axial  station  .  The  locations  were 
chosen  so  as  to  take  measurements  at  the  nozzle  exhaust  centerline,  at  the 
peak  velocity  location  and  at  a  radial  location  where  the  local  mean  velocity 


*  Note  that  there  was  a  year  lapse  in  time  between  the  conic  nozzle 

measurements  (performed  in  August  1976),  and  the  coannular  plug  nozzle 
measurements  (performed  in  August  1977).  During  that  interval,  improve¬ 
ments  in  the  data  acquisition  and  process  were  made.  Data  acquisition 
and  reduction  were  improved  by  a  factor  of  90. 

**  Note  that  at  each  measurement  station  (x,r),  several  axial  and  separations 
of  the  two  LVs  are  required  in  order  to  form  the  manifold  of  cross¬ 
correlations  depicted  in  Figure  10-5b.  In  general,  the  separation  distances 
used  were  e/Deq  =  0.04,  0.22,  0.39,  0.57,  0.74,  0.92,  1.09,  1.27.  Thus,  for 
each  location,  approximately  seven  cross-correlation  measurements  were  made. 
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Coannular  Plug  Nozzle  and  Test  Conditions  for  LV  Measured  Cross  Correlations. 


Table  10-2.  Summary  of  Turbulent  Structure  Properties  Measured  on  a 
Coannular  Plug  Nozzle. 


a.  Test 

-pi  = 
T 

Point  1  (P 

850°  R,  V* 

1 

£  =  1.8 

=  1250 

,  T°  =  1650° 
fps,  Vj>ix  = 

R,  V?  = 

J 

1602  fps) 

1758  fps, 

p£  =  1.8, 

X/Deq  r/Deq 

U,  fps 

U' ,  fps 

Vc,  fps 

Lx ,  in . 

vc/v,?ix 

C  J 

Lx/Deq 

1.81 

0.460 

1760.0 

96.84 

1587.6 

0.0495 

0.991 

0.0089 

2.54 

0.172 

1225.2 

63.71 

1227.1 

0.0288 

0.766 

0.00522 

0.241 

1370.9 

116.53 

1297.6 

0.1765 

0.810 

0.03199 

0.378 

1758.9 

91.46 

1754.2 

0.1715 

1.095 

0.03107 

4.71 

0 

1098.3 

121.9 

901.9 

0.2098 

0.563 

0.038 

0.329 

1606.3 

168.66 

1510.7 

0.3512 

0.943 

0.0636 

0.546 

1082.3 

284.6 

858.7 

0.5817 

0.536 

0.1054 

8.15 

0 

1228.0 

136.3 

1115.0 

1.0089 

0.696 

0.1828 

0.409 

1117.8 

271.6 

842.6 

1.1728 

0.526 

0.2125 

0.567 

879.5 

226.03 

765.8 

1.5775 

0.478 

0.2822 

9.97 

0.181 

1145.7 

221.1 

954.8 

1.3886 

0.596 

0.2516 

0.444 

1025.8 

245.2 

746.5 

1.4421 

0.466 

0.2613 

0.631 

816.0 

246.4 

730.5 

1.7545 

0.456 

0.3179 

b.  Test  Point  2  (P«  =  2.7,  T°  =  1650°  R,  V°  =  2216  fps,  p*  =  1.8, 
T*  =  850°  R,  V*  =  1250  fps,  Vj1X  =  1981  fps) 


X/Deq 

r/Deq 

U,  fps 

U' ,  fps 

Vc,  fps 

Lx,  in. 

vc/vfx 

Lx/Deq 

2.54 

0 

608.4 

119.8 

1010.0 

0.212 

0.510 

0.03845 

0.529 

1419.6 

359.0 

1010.0 

0.128 

0.510 

0.02313 

0.572 

981.9 

387.0 

1021.0 

0.283 

0.515 

0.05135 

4.71 

0 

1038.8 

133.9 

1295.0 

0.294 

0.654 

0.05334 

0.399 

2024.1 

151.8 

2009. 0 

0.272 

1.014 

0.04924 

0.546 

1704.2 

313.6 

1088.0 

0.687 

0.549 

0.1245 

8.15 

0 

1255.8 

251.2 

1338.0 

0.4649 

0.675 

0.08423 

0.314 

1456.2 

262.1 

1403.0 

0.9774 

0.708 

0.1771 

0.634 

856.11 

290.  3 

1887.0 

1.5105 

0.445 

0.2737 

9.97 

0 

1265.3 

270.8 

1271.0 

1.2423 

0.642 

0.2251 
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was  approximately  70%  of  the  peak  velocity.  Figure  10-10  illustrates  the 
typical  mean  velocity  profiles  for  test  point  1.  Illustrated  in  this  figure 
is  the  inverted  velocity  distribution  (high  velocity  on  the  outside  stream 
and  lower  velocity  on  the  inside  stream)  typical  of  this  type  of  test  nozzle. 
Also  shown  on  this  figure  are  the  approximate  locations  for  the  cross¬ 
correlation  measurements. 


Choosing  the  proper  normalization  or  similitude  parameters  for  the 
coannular  plug  nozzle  is  not  as  easily  done  as  for  the  conic  nozzle.  For  the 
coannular  plug  nozzle,  there  are  two  flow  streams  issuing  through  annular 
gaps.  For  the  turbulent  length  scale,  the  outer  stream  annulus  height,  h  , 
was  chosen  as  the  length  scale  characteristic  dimension.  For  the  convection 
velocity,  Vc,  normalization,  the  ideal  outer  stream  velocity  V°  was  used  for 
locations  X/Deq  <  5,  r/Deq  =£  0;  at  r/Deq  ~  0,  the  inner  stream  velocity,  Vj , 
was  the  normalizing  velocity  used  at  X/Deq  <  5;  and  the  specific  thrust, 
ymix,  was  the  normalizing  velocity  used  when  X/Deq  <  5.  Figures  10-11  and 
10-12  summarize  the  in-jet,  cross-correlation  turbulence  structure  measurements 
taken  to  date. 


Figure  10-11  illustrates  the  axial  variation  of  the  LV  measured  turbulent 

length  scale  for  the  coannular  plug  nozzle  at  three  streamlike  locations: 

region  of  peak  velocity,  70%  of  Local  U  ,  and  along  the  centerline  of  the 

max 

nozzle. 

The  physical  picture  which  emerges  from  these  measurements  is  the  length 
scale  of  the  convecting  eddies  are  smallest  along  the  centerline,  largest  at 
the  70%  Umax  stream  line  and  intermediate  in  length  along  the  Umax  stream 
line. 


Figures  10-12a  and  b  illustrate  the  convective  speeds  for  the  two  coannular 

plug  nozzle  test  conditions.  The  physical  picture  illustrated  by  these 

results  is  not  clear.  In  general,  it  might  be  postulated  that  the  larger 

eddies  moving  at  the  70%  Umax  stream  lines  are  convecting  at  the  slowest 

CH 

speed  (Vc  -  0.4  0.66  vV  ,  while  the  smaller  eddies  moving  along  the  center- 

line  and  at  Umax  are  moving  along  at  a  faster  rate  relative  to  the  chosen 
characteristic  velocity  (Vc  -  0.4  -+■  1.0  V^H) ,  with  the  higher  convection 
speeds  occurring  closer  to  the  nozzle  (X/Deq  <  5).  One  conclusion  which  can 
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Figure  10-10.  LV  Measured  Mean  Velocity  Flow  Field  for  Test  Point  1  In-Jet  Cross  Correlations 


Outer  Stream 


ial  Variation  of  LV  Measured  Turbulent  Length  Scale  for  a  Coannular  Plug  Nozzle. 
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(Open  Symbols:  Point  1;  Solid  Symbols:  Point  2) 


Figure  10-12.  Variation  of  LV  Measured  Convection  Velocity  for  a  Coannular 
Plug  Nozzle  at  High  Velocity  and  Temperature  Conditions. 


be  drawn,  by  comparing  the  conic  nozzle  data  shown  on  Figure  10-9b  to  that  of 
the  coannular  plug  nozzle  with  inverted  velocity  profiles,  is  that  the  convec¬ 
tion  speeds  downstream  of  the  nozzle  (X/Deq  >  5)  are  higher  for  the  conic 
nozzle  than  for  the  coannular  plug  nozzle.  Acoustically,  this  implies  that 
the  low  frequency  jet  noise  has  equivalently  less  convection  amplification 
and  fluid  shielding  (both  phenomena  are  strong  functions  of  convection 
velocity) . 

The  relative  turbulent  length  scales  between  the  conic  nozzle  and  the 
coannular  plug  nozzle  are  roughly  the  same  for  equivalent  X/D  (compare  tables 
of  Lx  in  Tables  10-1  and  10-2).  Smaller  eddies  are  close  to  the  nozzle  and 
larger  eddies  are  downstream  (high  frequency  noise  close  to  nozzle  exit  and 
lower  frequency  noise  downstream). 

10.4  TWO-POINT,  SPACE-TIME,  IN-JET  TO  FAR-FIELD ,  CROSS-CORRELATION  MEASUREMENTS 

During  the  same  test  period  in  which  the  two-point,  space-time,  in-jet/ 
in-jet  measurements  were  performed  on  the  coannular  plug  nozzle,  a  series 
of  in-jet  velocity  to  far-field  acoustic  cross-correlations  was  also  per¬ 
formed  with  the  LV  system.  The  primary  purposes  of  this  series  of  experiments 
were  to: 

1.  Identify  any  strongly  correlated  regions  of  the  jet  in  a  qualitative 
manner . 

2.  Evaluate  the  use  of  this  technique  as  a  potential  diagnostic  tool. 

The  general  test  arrangement  for  the  laser  velocimeter  was  the  same  as 
used  for  the  in-jet/in-jet  measurements  described  above  (and  described  in 
Section  9.0),  but  for  these  measurements  only  one  laser  was  used.  The 
microphone  used  during  the  experiments  was  usually  the  microphone  at  20°  to 
the  jet  axis*.  A  sketch  of  the  test  arrangement  in  the  Anechoic  Facility  is 
shown  on  Figure  10-13. 

*Note  that  to  perform  the  in-jet  to  far-field  acoustic  correlations  precisely, 
the  velocity  vector  in  the  direction  of  the  microphone  (as  described  in 
Section  6.0)  is  required.  For  the  current  measurements,  the  LV  control  volume 
was  set  to  measure  only  the  axial  velocity  component.  For  the  purpose  of 
these  experiments,  this  alignment  variation  with  respect  to  0^  =  20  was  not 
considered  serious  and  the  results  obtained  are  considered  to  be  representative. 


The  test  conditions  for  performing  the  in- jet  turbulent  velocity  to 
far-field  acoustic  cross-correlation  measurements  with  the  laser  velocimeter 
were  the  same  test  conditions  as  set  for  the  in-jet/in-jet  velocity  cross¬ 
correlation  measurements  described  above  (Section  10-3,  see  Figure  10-9): 
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That  is,  the  test  conditions  were  for  a  high  velocity  and  high  temperature 
coannular  plug  nozzle  with  inverted  velocity  and  temperature  flow  streams.  A 
summary  of  some  of  the  major  results  is  contained  in  Table  10-3. 

Contained  in  Table  10-3  are  the  maximum  values  of  the  normalizer  total 
cross  correlation  coefficient,  Rup*;  the  shear  portion,  Fshear '>  and  the  self" 
portion,  Rself.  of  the  total  correlation  function;  the  peak  frequency,  fp,  of 
the  power  spectrum  of  R-p,  and  the  maximum  of  the  second  derivative  of  R-p 
(near  ;  =  i0).  This  laser  ve  1  nr  inieter  measured  data  are  for  axial  and  radial 
point  locations  that  were  nearly  the  same  as  were  taken  for  the  turbulent 
structure  measurements  discussed  in  Section  10.3  above.  Figure  10-14  illus¬ 
trates  the  measurement  locations  on  the  laser  velocimeter  measured  mean 
velocity  flow  profiles  for  test  point  2. 

An  example  of  a  series  of  the  1,V  to  far-field  acoustic  cross-correlation 
measurements  is  shown  in  Figures  10-15a  through  g.  The  data  shown  in  Figure 
10-15  are  for  test  point  1  at  X/Deq  =  9.97,  r/Deq  =  0  (the  Appendix  shows  the 
radial  distribution  of  cross-correlations  also  found  for  this  case).  Figures 
10-15a  and  b  show  the  cross-correlation  functions  for  the  shear  and  the 
self-noise  portions  of  the  total  cross-correlation  function.  Figures  10-15c 
and  d  shown  the  total  cross— correlat ion  function  and  the  second  derivative  of 
the  total  cross-correlation  functions.  Figures  10-1 5e  and  f  are  the  laser 
velocimeter  measured  power  spectra  of  the  total  cross-correlation  function  and 
the  second  derivative  of  the  total  cross-correlation  function. 
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The  result;  of  Figure  1 0- 1  5  illustrate  features  which  exemplify  most  of 
the  measurements  taken  for  this  test  series.  The  values  of  the  measured  in¬ 
jet  velocity  to  far-field  acoustic  normalized  cross-correlations  are  of 
significant  magnitude  (for  this  illustrative  case  R-pot,-,]  =  0.261)  and  the 
actual  measured  correlation  function  shape  is  precise  enough  to  ensure  that 
"real"  regions  of  strong  correlati  ins  are  being  measured  between  the  laser 
and  the  acoustic  microphone.  Figures  10-1 5a,  h,  c  also  show  that  the  shear 
portion  (that  is,  the  mean/ turbulence  velocity  interaction  field)  dominates 
the  measured  total  cross-correlation  function.  The  preeminence  of  the  shear 
cross-correlation  was  found  to  he  true  for  100'  of  the  measurements  performed. 
Figure  10-15d  shows  the  second  derivative  measurements.  This  figure  illus¬ 
trates  that  the  new  statistical  estimator  procedures  developed  (see  Sections 
7.0  and  8.0)  together  with  the  improvements  o!  the  laser  processor,  were  success- 
fil*.  The  LV  measured  spectra  of  the  total  normalized  cross-correlation  function 
and  its  second  derivative  are  shown  in  hi  gates  1  < )  - 1 5  e  and  f.  These  spectra 
results  show  that  for  the  current  measurement  the  dominant  frequency  is  in 
the  lower  frequency  end  of  the  spectrum  and  is  fairly  clear  in  the  region  of 
importance.  For  RTotal  the  spectrum  is  quite  good  up  to  3  kHz;  thereafter, 
the  results  are  more  than  25  dll  down  from  the  peak.  For  the  spectrum  of  the 
second  derivative,  after  15  dB  down  from  the  peak ,  the  results  indicate 
further  improvement  could  be  made  to  define  the  higher  frequency  region. 

Figures  10-16a,  b  and  c  summarize  some  of  the  additional  results  obtained 
from  this  series  of  measurement! .  Figure  I0-I6a  shows  the  axial  variation  of 
the  maximum  value  of  the  total  in— jet  velocity  to  far-field  acoustic  pressure 
®jet  =  ^0  )  f°r  coannular  plug  nozzle  for  the  two  test  conditions 
tested;  the  results  are  for  the  centerline  measurements  only  .  The  values 
of  RTotal  show  a  range  between  0.029  to  0.261,  depending  on  axial  location; 
with  the  maximum  occurring  between  8-10  IJeq .  For  test  point  2  (the  higher 
velocity  coannular  plug  nozzle  test  case).  Figure  10-16b  illustrates  that 
Rshear  generally  20  times  greater  than  R.self  for  all  of  the  measurements. 


These  new  results  are  an  up-date  to  previous  results  reported  in  Vols  I  and 
II  of  Reference  31. 

ick 

Peak  values  of  were  found  along  the  centerline.  See  Table  10-  1  for 

values  at  other  radial  locations.  Also  see  the  Appendix. 
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Figure  10-16.  Laser  Velocimeter  In-Jet  to  Far-Field 
Measurements  for  a  High  Velocity  and 
Temperature  Coannular  Plug  Nozzl®. 
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Figure  10-16c  illustrates  the  axial  variation  of  the  peak  frequency  for  both 
nozzle  conditions  tested.  These  results  show  that  the  measured  peak  frequency 
of  the  Four ier-Transform  of  K'|otal  decreases  with  increasing  axial  location 
(higher  frequency  close  to  the  nozzle,  lower  frequency  further  downstream). 

The  rate  of  decrease  in  frequency  is  somewhat  slower  tnan  what  would  lie 
expected  from  a  conic  nozzle.  For  illustrative  purposes,  it  was  found  that: 


i  I)e<i  _n  s  7 

—  -  ( ri .  28  X/Deq) 

vmix 

i 

is  somewhat  representative  of  the  measured  findings. 

As  an  example  of  LV  measured  acoustic  intensity  per  unit  volume.  Section 
6.0  showed  that  the  quantity  to  be  evaluated  was: 


1/  AV 
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A  measure  of  this  quality*  is  illustrated  using  results  from  test  case 
2  at  X/Deq  =  9.97,  r/Deq  =  0.  The  maximum  value  of  •  >  “^u^p  ’  /  fit  2  was  found  to 
be  3.2  x  1015  pbar-f t2/sec^ .  The  distance  from  the  measurement  station  to 
the  acoustic  microphone  was  24.27  feet.  Thus,  the  computed  acoustic  intensity 
level  at  fp  ~  500  Hz  is  found  to  be  -  110  dB  at  0j  =  20°.  Measurements  of 
acoustic  intensity  for  this  nozzle  at  the  same  VjUX  (Reference  52)  was  found 
to  be  114.8  dB  at  f  f / 3  ~  500  Hz.  The  computed  value  of  110  dB  from  the  LV 
measurements  is  considered  to  be  quite  representative  since,  in  fact,  it  is 
expected  that  the  500  Hz  spectrum  is  composed  of  noise  elements  from  several 
jet  volume  sources. 

The  value  used  to  approximate  S^u^p ' /9x ^  I  was  (-f,e  vapue  Gf  3  ^u^p'/Tit^ 

at  1 1  io  t  which  Ruzp,  was  measured  as  a  maximum.  This  quantity  is  called 
D4;  in  Table  10-3  and  is  also  shown  on  the  cross-correlation  figures.  The 
results  have  shown  that  the  t  for  which  Ru  P  is  a  maximum  always  occurred 
slightly  prior  to  the  calculated  t0.  The  reason  the  computed  t0  value 
itself  was  not  used  was  because  the  real  retarded  time  for  these  test  cases 
are  known  to  be  less  than  the  simple  t0  =  r/ac.  The  reason  this  is  true  is 
that  part  of  the  sound  transmission  path  is  through  the  heated  portion  of  the 
jet.  During  that  period  the  sound  is  transmitted  faster  (due  to  the  jet's 
higher  temperature)  than  the  ambient  speed  of  sound.  Thus  the  actual 
retarded  time  should  be  expected  to  be  somewhat  less  than  the  r/a0 
calculation . 


1  1  B 


The  results  illustrated  above  have  shown  that  the  laser  velocimeter 
techniques  developed  over  the  past  several  years  can  be  useful  as  a  diagnostic 
tool  for  better  understanding  realistic  jet  nozzle  turbulent  and  acoustic 
features  at  typical  engine  velocity  and  temperature  conditions.  The  results 
have  illustrated  that  important  aerodynamic  source  term  results,  such  as 
turbulent  length  scale,  convection  velocity,  turbulent  and  mean  velocity,  can 
be  measured  for  complex  nozzle  flow  fields  and  that  cross-correlation 
measurements  between  the  point  velocity  field  and  the  far-field  acoustic 
radiation  field  have  shown  strong  regions  of  correlation  heretofore  not 
measured.  These  results  are  considered  to  be  quite  encouraging  for  future 
investigations  in  nozzle  noise  and  aerostructure  characteristics  of  high 
velocity  and  temperature  jets. 

Some  additional  remarks  regarding  the  cross-correlation  function  should 
be  noted.  As  shown  in  Figure  10-15c  the  measured  normalized  cross-correlation 
function  is  somewhat  symmetrical  with  a  negative  peak.  This  shape  was  not 
always  the  same  for  all  measurements.  For  certain  locations,  the  cross¬ 
correlation  coefficient  was  symmetric  with  a  positive  peak.  An  example  of 
this  type  of  correlation  function  is  given  in  the  Appendix  (X/Deq  =  9.97, 
r/Deq  =  0.246).  The  reason  for  such  a  contrasting  difference  between  the 
two  cross-correlation  functions  is  expected  to  be  related  to  some  physical 
flow  phenomenon  which  needs  further  evaluation  and  explanation. 
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11.0  CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


11.1  CONCLUDING  REMARKS 

The  studies  and  discussions  contained  in  this  report  include  a  number  of 
unique  and  modern  concepts  for  laser  velocimeter  signal  processing  and  data 
handling  and  reduction  associated  with  performing  real-time  velocity  and 
cross-correlation  measurements  in  realistic  jet  exhaust  nozzle  flows.  The 
introduction  of  the  filter  bank  approach  for  LV  signal  processing  and  the 
use  of  quantized  product  factoring  in  the  LV  data  handling  cycle  have  enabled 
a  time  improvement  of  a  factor  of  90  in  LV  data  acquisition  in  the  most 
severe  type  of  velocity  and  temperature  conditions. 

A  series  of  two-point  space-time  velocity  measurements  for  determining 
a  jet's  turbulent  length  scale  and  convection  speed  was  successfully  per¬ 
formed  for  a  baseline  conic  nozzle  and  for  an  inverted  velocity  and  tempera¬ 
ture  flow  coannular  plug  nozzle.  Experiments  were  performed  in  jet  exhausts 
of  2200  fps  and  1650°  R.  Additionally,  a  series  of  in-jet  velocity  to  far- 
field  acoustic  cross-correlation  measurements  was  performed  on  a  coannular 
plug  nozzle  at  sonic  and  supersonic  heated  flow  conditions.  Measurements  of 
this  type  have  not  been  previously  performed.  The  results  of  the  studies 
performed  to  date  indicate  that  the  LV  technology  is  at  hand  to  begin  sys¬ 
tematic  aerodynamic  and  acoustic  noise  source  type  experiments  for  high  tem¬ 
perature  subsonic  and  supersonic  exhaust  jets. 

Some  of  the  results  obtained  from  the  recent  cross-correlation  measure¬ 
ments  are: 

•  For  a  simple  conic  nozzle,  there  exists  a  radial  and  axial  vari¬ 
ation  of  turbulent  length  scale  and  convection  speed  throughout  the 
jet. 

•  For  the  conic  nozzle, it  was  found  that  not  a  great  difference  be¬ 
tween  ambient  and  heated  jet  flow  existed  for  length  scale  and 
convection  speed. 

•  For  the  inverted  flow,  coannular  plug  nozzle  tested,  an  axial  and 
radial  variation  of  length  scale  and  convection  speed  also  was 
found.  The  physical  picture  that  emerges  is  that  the  smaller 
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turbulent  eddies  of  the  exhaust  plume  are  convecting  along  the 
center  line,  the  larger  eddies  are  at  a  70%  Umax  stream-line  on 
the  outside  of  the  jet,  and  intermediate  size  eddies  are  convecting 
along  of  the  Umax  stream  line. 


The  convection  velocities  for  the  coannular  plug  nozzle  were  found 
to  be  highest  close  to  the  nozzle  and  to  gradually  decay  in  speed 
with  increasing  axial  location.  The  centerline  and  Umax  streamline 
eddies  convected  the  fastest,  with  the  slower  convecting  eddies 
located  at  the  70%  Umax  streamline.  Typical  convection  speeds 
ranged  as  Vc  =  0.4  —  1.0  Vcharacteristics • 

Regions  of  strong  correlation  were  measured  on  the  coannular  plug 
nozzle  where  the  LV  velocity  signal  was  cross-correlated  with  an 
acoustic  microphone.  Normalized  cross-correlations  of  0.26  were 
obtained. 

For  all  tests  performed,  the  shear  turbulence  aeroacoustic  correla¬ 
tion  function  dominated.  The  shear  cross-correlation  coefficient 
was  at  least  a  factor  of  20  greater  than  the  self-turbulence  aero¬ 
acoustic  cross-correlation  function. 

The  peak  frequence  of  the  cross-correlation  function  varied  as 


f  Deq 


V' 


.mix 


=  (5.28  X/Deq) 


-0.57 


Thus,  the  classic  notion  of  high  frequency  noise  generation  occurring 
close  to  the  coannular  plug  nozzle  and  low  frequency  noise  occurring  down¬ 
stream  was  verified. 


11.2  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  The  development  of  the  laser  velocimeter  as  a  viable,  economical 

jet  noise  diagnostic  tool  has  been  demonstrated.  It  is  recommended 
that  a  series  of  systematic  and  detailed  measurements  be  performed 
to  better  define  the  turbulent  structure  of  jets,  convection  speed 
and  scale  of  turbulence.  The  initial  nozzles  recommended  for  study 
are:  1)  the  conic  nozzle,  2)  coaxial  coplanar,  3)  coannular  plug 
nozzle.  Experiments  should  be  performed  at  subsonic  and  supersonic 
conditions.  In  addition,  in-jet  velocity  to  far-field  acoustic 
pressure  cross-correlation  measurements  on  the  above  three  nozzles 
should  be  performed  to  quantify  the  regions  of  strong  correlation 


Application  of  some  of  the  more  modern  theoretical  acoustic  model¬ 
ing  techniques  toward  additional  definition  and  clarification  of 
the  aerodynamic  source  terms  should  be  performed. 

Although  great  strides  have  been  made  in  the  improvements  of  the  LV 
processor  and  LV  head,  additional  advancements  will  be  forthcoming 
with  the  applications  in  large  scale  integrated  circuits  (micro¬ 
processors)  .  Development  as  a  Fast  Fourier  transform  type  LV 
processor  with  a  microprocessor  approach  is  recommended. 

Expansion  of  the  LV  arrangement  to  include  a  multivector  output 
arrangement  is  desirable. 


APPENDIX  -  LASER  VELOCIMETER  IN-JET  VELOCITY  TO  FAR-FIELD  ACOUSTIC  CROSS- 
CORRELAT ION  MEASUREMENTS  FOR  A  COANNULAR  PLUG  NOZZLE 


•  Measurements  Taken  at:  X/Deq  =  9.97;  r/Deq  =  0.246, 

•  D_„  =  5.519  inches 

eq 

•  Test  Conditions: 

-  Outer  Stream:  Pr  =  2.7,  T^  =  1650°  R,  Vj  =  2216  fps 

-  Inner  Stream:  Pr  =1.8,  TT  =  850°  R,  Vj  =  1250  fps 


0.661 
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RTotal  RSelf  RShear 


RTotal/^T2  RTotal 


S(32RTotal/3t2>  S(!W> 


Total  self  KShear 


Total'81'*  S<RTotal> 


Total  self  Shear 
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